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1.0 INTRODUCTION 


Volume 2 of the final report is used to outline four computer 
programs which were developed during the grant period. All were coded in 
FORTRAN IV Language for usage on the University ofVhshington CDC 6400 
digital computer. 

The discussions here are primarily user input instructions. For two 
of the programs, ROTOR and AIRFOIL, the background models are developed 
in some detail in Volume 1 oi the report. Program TONE is developed herein 
and SDATA is outlined here with reference to the literature where the 
detailed development requires a book to .describe. The function of each 
program is discussed b€:low, 

(1) Program AZilFOIL; This program computes the spectrum of 
radiated sound from a single airfoil immersed in a laminar 
flow field. The m. chanism is force fluctuations which are 
related to the unsteady wake momentum. Input required is the 
aerodynamic properties of the turbulent wake , 

(2) Program ROTOR; The purpose of program ROTOR is to provide an 
extension of the single airfoil in AIRFOIL to a rotating 
frame. This is a model for sound generation in subsoni,c 
rotors. The only broad band mechanism currently in the pro- 
gram is the airfoil wake mechanism. Program ROTOR also 
computes tone sound generation due to the steady state forces 
on the blades (the Gutin result) . Input requirements arc 


similar to those for AIRFOIL, 
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(3) Program TONE; A moving source analysis Is used to generate 
a time series for an array of forces moving In a circular 
path as would be the case for a rotor. This analysis has 
great flexibility In that arbitrary blade spaclngs, blade 
forces and source motions can be easily accommodated. Broad 
band effects have been Included using autoregression methods. 

The resultant time series are Fourier transformed using a 
Fast Fourier Transform to present the sound radiation in 

the more normal spectral form. With some study, a user can 
make this program a very valuable tool in sound generation 
Investigations . 

(4) Program SDATA; This program is a standard time series analysis 
package. It will read in two discrete time series and form 
auto and cross covariances and normalize these to form 
correlations. The program will then transform the covariances 
to yield auto and cross power spectra by means of a Fourier 
transformation. The spectral data are presented In terms of 
amplitude and phase as well as the coherency spectrum. 

A final note on these programs is in order. These were developed as 
tools to be used in the work performed under the grant. There was no 
intention to provide the programs as complete packages for general usage. 
Tlie content is general enough, however, that with some effort a user should 
be able to employ the programs to advantage. He is cautioned to understand 
the analytical models before he becomes too involved in turning around 
decks. No program is a substitute for this understanding. Tools must be 
skillfully used. 
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To the end of providing a usable tool the programs are coded in 
variable names which read very like those of the physical variables they 
represent. For example the density will be RHO and the longitudinal 
turbulence covariance will be UlUl, etc. This should facilitate the 
reading of the program listings. 

Each of the programs is discussed in the following sections. 


2.0 PROGRAM AIRFCIL 


The fortran yrogram AIRFOIL is used to compute the sound radiated 
from a single airfoil due to those mechanisms related to the airfoil wake 
turbulence. A development of the model was presented in detail iii Volume 
1 of this report. Through the model, the sound radiation is related to 
the turbulent structure in the wake of the airfoil. 

The equation which forms the basic result of the model and that pro- 
grammed here is 

2 2 

G (x’,(u) =5] 16 S L^(u u , m) L, (u^u ,(o) L^(u^u ,(o) 

P i=l j=l J *• J J 

10 A A u u (y„,to) dy„ 

with: A, = p U, sin i|j/CA'ira^ lx*") 

1 o i o 

A- = p U, cos t|//(4Tia^ lit I) . 
z o 1 o 

L .menclature used here is the same as that of Volume 1. 

The integration noted above is performed using a Lagrangian method. 
This subroutine is part of the University of V^ashington Computer Library. 
Any simple numerical integration scheme may be substituted for this sub- 
routine. In the program the subroutine is labeled LAGRAN. 

2.1 Program Input 

AIRFOIL is a simple computational program. It consists of three sub 
programs. The main program is AIRFOIL in which the actual computations 
are made. The next is subroutine DREAD in which data are read into and 


conditioned to form the variables used Iji the actual computations. Finally 
DWE.ITE in which the results of the computations are written out. 


Input to the program is accomplished using IBM cards. In the follow“ 
ing we ^<111 list each card type, the variable names and their descriptions. 

CMD 1, FOEMAT(7F10.2) 


PATH 

atmosphe: c static pressure 

(atm) 

TATM 

" " teinperature 

C°c) 

SPAN 

airfoil ength 

(cm) 

XLAM 

length in the stream direction (this is 
a dead variable in the current program so 
any input is okay 

(cm) 

DEL 

wake thickness 

(cm) 

RAD 

distance from airfoil to sound observer 

(m) 

CARD 2, FORMAT (215) 


NFF 

number of frequency points where data will 
be input 


NPY 

number of positions in the wake where data 
will be given 


CARD 3, FORMAT (7F10. 2) 


Y(I) 

I = 1,NPY 

positions in the wake measured relative to 
the wake centerline 

(cm) 


Note that if more than 7 entries are needed, additional cards will be 
required for this variable. 

CARD 4, FORMAT (7F10. 2) 


FREQ (I) 

frequencies at which the following 
correlation lengths are to correspond 

(Hz) 

XL(I) 

correlation length in stream direction for 
the longitudinal turbulent component 

(cm) 
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YL(I) 

ZL(I) 


correlation length across the wake for 
the longituf^lnal turbulence component 

correlation length In the span direction 
for the longitudinal turbulence component 


CARD 5, FORMAT (7F10. 2) 


FREQCD 

same as for CARD 4 





B 

xr(i) 

same as for CARD 4, but 
turbulence component 

for 

the 

transverse 

(cm) 

ii 

r ! 
.It 

YT(I) 

same as for CARD 4, but 
turbulence component 

for 

the 

transverse 

(cm) 

■\ i 
-.1 
" ! 
• i 

i\ 
- 1 

ZT(I) 

same as for CARD 4, but 
turbulence component 

for 

the 

transverse 

(cm) 

; L 

^^4 


Notice that CARD’s 4 and 5 are read inside a DO loop. The statements 
are 

DO 16 I = 1,KPF 

READC5 »3) FREQ(I) ,XL(I) ,YL(I) ,ZL(I) 

16 RFAD(5,3) FREQ(I) ,XT(I) ,YT(I) ,ZT(I) 

Certainly more than just the two cards will be required here since more 
than one spectral point will generally be required. For this outline, we 
will continue with CARD 6 keeping in mind that CARD 6 riifers to a new input 
type of data rather than the actual number of the card in the input file. 


CARD 6, FORMAT (7F10. 2) 

XI (I) ratio of wake velocity at Y(I) to the velocity 

1 = 1,NPY in the free stream UMEAN 


CARD 7, FORMAT (7F10. 2) 

DDBl(I) ratio of spectral level at FREQ(I) to the RIIS 

1 - 1,NPF level; this is for the longitudinal component 

of turbulence 
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DDBl can be interpreted from the equation for its usage which is 
uj^(f) = * 10. **C“DDB1(I)/10.) 

u^(f) is the actual spectral level sought, and 

u^ is the RMS level for the longitudinal component of turbulence. 

CARD 8, FORMAT (7F10. 2) 

DDB2(I) the ratio of the spectral level at FREQ(I) 

I = IjEPF to the RMS level for the shear component 

of turbulence (dB) 

see the note for DDBl 

CARD 9, FORMAT (7F10. 2) 

DDB3(I) the ratio of the spectral level at FREQ (I) 

to the RMS level for the transverse component 
of turbulence (dB) 


CARD 10, FORMAT C7F10. 2) 


UllCD 

I = Ij^R’Y 


spatial distribution of the RMS level of 
longitudinal turbulence across the wake, 
normalized with UMEAN 


CARD 11, FORMAT C7F10. 2) 

spatial distribution of the RMS level of the 
shear stress across the wake, normalized 
with UMEAll 


CARD 12, FORMAT C7F10. 2) 

U22(I) spatial distribution of the RMS level of the 

transverse turbulence across the wake, normal- 
ized with UMEAN 

This completes the inp r_ for the program. 


UlKD 
I = 1,JIFY 
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Program output consists of a listing of the input variables and of 
the spectral distribution of the radiated sound pressure at each radia- 
tion angle 0° through 180° in 10° increments. The radiated sound is 
presented in dB by normalizing the pressure i>rith a reference 20pN/in^. 


ririririooriornr'innnnnnnoononn n n n o 


2;i ProMia likiiiia 

HKOGRAi'-’t Al'<F<jlL i iNHUT ».uuTPUT tTAPhS-lNPuT* 1 aPE6=^UTPUT ) 

THib program COMPUTtb The GRtOTRUM OF RAUIaTlu 50UND FROM 
■ CORRcUAXiM : AN bpELTRAE VATA ^ POINT VtLOOi I Y 

CORpt-tATldRS TN TKE ;>mAKt: oK a AiRFoiL 

DIMENSION XINTGdl) »XlNTl ( 11) »XINT( XI ) 

D i M £N b I ON . F R Ew (20 } ,XL(^:0) TYLl2v>|*iL('2U)».Uiul {11.2 V V »X1 ( 11) 
10102(11.20) »v2U2(il»*:U) .Y (11) .Xi-idl) .XunEG ( 2u ) . 

2RbD(2V) ,:3 Rl12v) »v^cLl2v) »vCi'(2v^) »VAk( 3) *X1 ( 2o ) 

DATA PI/3.1AiS9/ ; 

READ IN DATA . 

P-ATM = STATIC pressure (MM hO) 

TATM = static temper a tORE (DEb C) 

Span. = airfoil length ( cm ) 

XLAM = STREAM INTEGRATION LENGTH (CiM) 

DEL = WARE THICRNESv (CM) 
uMeAN = FREE stream VELOCITY (h/SeC) 

RAP = DISTANCE TO OBSERVER (M) 

NPF - Nui^oER OF Frequency points in speltrum 
NPY = NUMdlr OF SpATlAL POINTS In WAKl 
Y = EOCATTONS OF THE POINTcj IN THE WAKE (CM) 

FREQ = FREQUENCY ARRAY (HE) 

XL = stream correlation For ulUl (CM) 

YL = CORRiiLATION ACKOSo I Ha »<nKii EOK ulul (CM) 

AL = SPAN CORRELATION FOK UlUl ( C.V| ) 

XT = stream CORRELATION FOR U2u2 (CM) 

YT = C0RR<-LATi0N ACROSS THE wAKt eOR U2o2 ( CM ) 
n = SPAN correlation FOK U2U2 (CH) 

XI = RATIO OF WAKE VELOCITY TO UMEAN ' 

DDBJ. = RATIO OF SRECTRAL LEVEL 70 KMS - UiOl (DB) 

DUbZ = RAllO OF SpECTRAL LE\/LL Tq RMb - ulv2 lUo) 

ppo3 = RaIIo OF optCTRAL LLVtiL lO kMo - u2U2 (Do) 

vJli = RMo value Of UlUl (PEKCENI) 

ul2 = RMS value of U1U2 (PERCENT) 

022 = RMG VALUE OF U2U2 (PERCENT) 

CAL L DR EAU .( A(J » RHO , P I » BM » XOM EG . V C L » V C 1 » G PAN » X L AM » D EL » X M , 

IK AD » ulul » ^102 » U2o2 » NPF * NPa * itP i >X » t » F Khw * KATl''l » ! A 1 M » AL » aT ) 

C COMPoTh RavIaTION PATTLRN 

COEFl= ( RH0/(4.->^PI*RAD) ) **2 
C0EF2=C0EF 1*SPAN*3M**2 
DO 2000 KK:=1,19 
UKK-KK 
C0DE1=U1CK 

THETA=-9U.+HJ**(UkK~1.) 

THET=THETA*PI/18V, 

ST=SIN(IH£T) 

. CT=:C0S(THET) 

STS-ST»*2 
.CTS=CT*»2 
5TCT=ST4CT . 


EEPKODUCIBILITY OF THE 
QEIGIlSrAL PAGE IS POOR 


10 


;C ' eOf^PUTE iPECTKUivr " . 

. ' Dd'^lOpO <=1»NPF ^ . ,■ '■ 

^C- 'ClJMPUTE Inteor^ over the ^mke volume . 

XTNTt<:)^0*U . , 

DO 30U il=i>3 
DU; 100 I=x>:MpY . 

VAR ( 1) =i>T.i>*UlUl { N. » I ) **2 *VcL( s ) *XL{ Is) 

^ yAl^tijT=Stcr^^^^^ ) /2.*UL(A)+P<Tt A) )/2. 

: VARt3i=CTo^U2U2TK>r)^^^ 

100 X.Ii^T6 ('I)-VAR( 11 ). 

^ wRllE^^ 

10 format ( lEl^. 2) 

call LAGR/^N( V >XINTG »(NPY »i >1 jAin IkKOR ) 

;30g XlivTCK) = AJ'*S+XlNT{K) 

pSDLK)=XlHT{K)*COFF3*XOMt;G**A 
PSD{^)=ABb(PSD(K.) 1 
WRTTEt6*ll) ANS 

11 FORMAT (4H ANS 1e2^*5) 

1000 oPu(n) = iO,*aLug1'^(PSD1 A)/4ou.e-12) 

gall DWRI 1L ( AO*RriO»Pl »BM»XUrtLb» vCL»VCT •^HAn»XLAi-i»DEL » 
IXcUKADtUlOi ,0iU2;»U2o2 ♦XiYV THE T A »PbD » GHL » HP Y i NHX » i'<PF » 
2FREG»CODEi»PArM*tATM»XLfXT) 

2000 CONTINUE 

or OP 
END 

bUbRuuT iNb DREAD I AU * KriU »H I » Bi^i » XUMtO i VGL » VCT »oPAH» XEAis » DEL jiXri » 
1KAD»U1U1 j«JlU2 »U2U2*iwPFtHKX»JMKif »X. Y » fREU »PATci » FA f »XL »XT ) 

C THlb SUBROUTINE READS ANU CONDITIONS INPUT DATA 

UIMENSION FR£Q(2G) .XL(2y) *YL(2w) Zw ) *ulul i 11 *2y ) .XI Ul ) . 
lJiU2(11.2y ) »U2U2< 1I»20) .r III) iXi'Kil) .AUnEGUo) » 

2 PsDtZu} .oPLlZU) »VGi-(2y) »VCI t2y) 
uIMlNSION DDBitZy )>DDfcS2(£y) .DDbJiiu) tUlK lU *olii( il> • 

1^22 (11 ) »DF1 ( 2U ) *DF2 ( 2D) .DF3( 2^ ) *Xl ( 2^ ) » Y T ( 2u ) .2 l( 2D ) 

DATA CR.Gu»GAM/2a8. .1.0V1.4/ 

READtS.l) PATM.TATM*DPANtXLAM.DEL*UMEAN.KAD 
1 FORMAT (7F10. 2) 

CMFT = ,0320 

PATM * 2li6.*PATM 

T atm* ( TATM+32 • ) *9 ^ /5 #+460 . 

SPAN*SPAN*CMFT 
. XtAH»XLAM*CMFT 
DEL»DEL*CMFT 
UMEAN T UrtEAN*CMFT*100. 

RAD » RAO*CMFT»lOo. 
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dscai:e=i*o 

rK0=RA.TM/ 

REAG(5iv2) NPF^NPY ' ■- 

2 ■fOREA:T:G3:i 5 V^ 

: ' READC-5*!) iyil)»T = X»'NPY) 

DO- 18 I=1»NPY 
18 Yil)=CWFT*Y(;i) 

. DG 16 :I = l»^iPF 
; READI5V31 

16 READ t 5* 3 ) :. FReq ( I ) »KT C I ) jYT( I ) *^-T t I ) 

.DG: .1Q; K=1*NP:F 
XOHEG(K'y=2*#I^ 

:FR:EGlKi=FR£GlX)/.DGCAL£ 

X6.yiEGt(Ci = 

GGtF-2.w/l ^:..54:i^l2v ; / 

xG(Kl=XEli^HcbEF; 

Y LI fc)^=Y L tx y.#GGEF- 
■y : ZLiNy =zGi;^y#^^ 

XTCR^=XT{rs)XCQEF . 

: YT|Ay=VT(N)*COEF 

ArtR)=ZTiiG)/XGO,E,F 
. vGL:{:(^i~XLiK):*YLlW 
10 VGXlfc)-Xiri<.):^YTiX}'K^^ 

READ( 5 »1) (XI ( I ) *I-1»NPY ) . 

. \D0 13 I^lrNPY 
13, XM( 1 ) =OMEAMXi ( l) /AO 

READ(5*1) (DDtil(Rl»K=l»NPF) 

: isEAOC5»l} (0UD2iP^i*A=l»MKF). 

K E A D C 5 * 1 ) ( D Dd 3 ( K. y * K = 1 ♦ NP F ) 

DO 12 K=1»NPF 

DF1{ Ay=lO,«-*(-DDBl < /2'J. ) 

DF2tR) =lO.F*( -DDB2 Fi<- ) /2i^ 

\Z DF3(X) = 10.**t-DDB^ f X}/2'J. ) 

READ(5»1) {.‘^liny»I = l»NPYl 

READ(5*1) .(UI21I ) ♦I=l»NPYj 

REAG(5»13 (U22 I I ) *:I = l»NPYl 

DO 15 i<=i»NPF 

DO 15 I=1,»NPY 

uun r )=DFUA)*U11 ( I) 

J1U2FK* n =DF2t K. )*01Gt I ) 

15 O202fK*Il=DF3{'<).*U22(n 

3 F0R|V.AT ( AFlu,2) 

RETURN' 

END : : . 

oUBR'-'UTTNE DWRI TE ( A'J » P I » bi*t » AuriEG »VCi->VCT »6 i-‘AN * XL Ah s DEl s X(V, j RaD » 

luiui si UlUZ »u2U2 »X* Y * THET a *PBp »bPi- ,ivH Y »iiPX > NPF » F^E'^ »C0DE1 » 
2PATRi*TATN*XL»XT] 


5SSi'?S°- 


L)I.-lt!‘{SIUN FKEu(2uj »XL(2uT,YLt2o),/lU(2v/}»Ului(li»2'J)*XMll)» 
loXU2(ll92'j)»02U2Cllt2j) sX(''i (XI) * AUi''itG t 2»J ) » 

2 PoD(2o} »oML(20) » VCL t 2*^ ) »VC T ( 2^ ) 

IF(C0DE1«GT,1. ) GO TO 21 

l/'*RlTE(6il ) 

1 HoRMAT ( IHi » 8^H bOUNL' RAUIaIIOM FROM a biNGLt AIkF'oIl AS RciLaTlj TO 
l!/.AiCL AEROOYNAivlIC pARaMcTlKS //) 

rtKiTE (6»2 ) PAfFi* I AT''i»bPAN*XLAn»iJEL»KAO 

2 FORi''iAT ( 1 7ri AYidlEiNT PKEbSuRL lU ,/L » 2'^h At'iblENT TcMPERATuKE 

llFi 0 . 2 // 2 on blMc-NSlOMo OF INTuoKaIiON /13M AlRrOiL oRaN 
2it20.5»X0X.17H stream ulRcCTiON aAKl THICnNLSS 

3iL20.S//iy., RADIOS TO Od^ERvLK ic2o,2///) 

rtkl TE ( 69 3 ) ( Y ( I ) 9 1 = 1 »NPY ) 

3 FoRjV,MTCbUn STlauY FLU/. NOMltlK i Tkc/vM jl R(_t„ 1 i OM 

13oX9j.6H A'^KUS^ TmE v'.AaE / * i o F X • A/ ) 

wKITE( 69A) tXM( I ) 9 1 = 1 »iMFY) 

A FoRi'.AT ( 5X® 1F10.A» 5 X 9 lOFlu.3 ) 

DO 20 <=l9NPF 

wRnE(6»6) FREQ(k) » v CL ( K ) > v 0 1 I j 

6 FORi'lAT ( irix 9 luM FREO'JtNC Y lL2u8b»lVri COiiK VOLu.-tt (LONG) iL^^e:;* 
119h CORK VOLOi'iE (TRAN) l£2-'0//) 

HRirE(69 7 ) ( Y ( I) , I = 1 9NHY ) 

7 FORMAT(:jAH LUNG I TUN I NAL ToRdULtNCL CUjV|RUNENT //i7h STREAM DiRuCTiO 
IN JUX^iuH ACRUbS 1 HE w aKL/ i^A » luF 1 '-i o 2 / ) 

rtRITL(69A) ( 0101 ( K 9 L ) 9 1 = 1 9 NPY ) 

rt^IiE(6»8) ( Y { I ) 9l = l»NPY ) 

8 FUKhA T { / / Jilri i RANSVcKoE TuKcJ^jlLNCc. CumRui'^EN i / / 1 V i ) .j Tk[_A>-i Lj Ikc'-T 1 

IduXsibri A'-KUbo TnE /2^A 9 i uF lo • 2 ) 

rtRlTE(6*A) (0202 ( tv9 I } 9 1 = 1 9NrY ) 

rtRlTE(699) (Y( I) , 1=1 »NPY ) 

9 Format ( //-i7n oNC-ar TurduLi-NOL COMPOiNONr //17n oTRcam uiRuCIiON 
liUX^ibFI AoKObb TnE viiAXE /2 JX 9 XuF lu . 2 / ) 

i«RnE(69A} (U1U2 I K 9 I ) 9 I = I iNPY ) 

20 CONTINUE 

21 CONTINUE 

wRlTE(69ll) THElA 

11 FORMAT ( irU »25h PREOlCTEO ^OuNo optiCTRoM //Ibti RAuIATiON ANOLt- 
UFi0.2//10ri FREuueNCY ioXfl/n kAOIAThl) bPL(ub) ) 

I 

DO 12 I=1*NPF 

12 WRITE(69i3) FREOmsSPLd) 

13 FORMAT! IFlOaZtlUXjlFlQ^E! 
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RETURN 

ENu) 
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3.0 PROGRAM ROTOR 

The purpose of developing program ROTOR was to provide a method for 
the prediction of sound generation from turbomachinery rotors. Aerodynamic 
operation was limited to the subsonic flow regime. Both tone noise and 
broad band effects were included. The tone noise part of the model was 
simply the steady state forces on the rotor blades. The broad band mech- 
anism was taken as the force fluctuations which were relatable to the 
wake turbulence of the individual airfoils of the rotor. The siathematical 
model used for the basis of the sound radiation was that presented by 
Ffowcs Williams and Hawkings'. 

To use the program one must supply the aerodynamic parameters for the 
particular design. Program ROTOR then computes the spectrum of sound 
radiated by the machine due to the mechanisms noted above. 

The program is divided into five subroutines. A brief discussion of 
the contents of these subroutines will be given to provide some familiari- 
zation with the program structure. The names of the subroutines are; 
DATAIN, WAKE, ORTHl, TONE, and 0RTH2., 

3.1 Subroutine OATAIN 

As the name implies, DATAIN is used to read and condition the input 
data. The data input can be listed with the required FORMAT and 
physical units. The following are the only inputs required to run the 
program. 

CARD 1, FORMAT (7F10. 2) 


RPM 

rotor speed 

(RPM) 

R1 

distance to observer 

CM) 


"Theory Relating to the Noise of Rotating Machinery," J. Sound if Vib., 
1, 1969. 
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B 

number of rotor blades 

C-) 

01 

blade relative velocity 

(M/SEC) 

XLAM 

blade stagger angle 

(DEGREES) 

RS 

source radius 

(M) 

FMAX 

maximum frequency in spectrum 

(iiz) 


CARD 2, l?omT(7F10.2) 



FB 

magnitude of force on each blade 

Cltg) 


UlUlB 

RMS turbulence level for the longitu- 
dinal component, normalized with U1 

(— ) 


U1U2B 

RMS shear turbulence, normalized x<ath U1 

C-) 


U2U1B 

RMS shear turbulence, normalized with U1 

(-_) 


D2U2B 

RMS turbulence level for the transverse 
component, normalized with U1 

C-) 

CARD 

3, FORMAT (7F10. 2) 



ZL 

blade span or other span length used to 
define the length of the source in the 
span direction 

(M) 


YL 

thickness of the airfoil wake 

CM) 

CARD 

4, FOR>tAT(lI5) 



NF 

number of frequency points to be used to 
input random data for wake related sound 
generation 

(-) 


In the remaining inputs a FORMAT (7F10 ,2) is used but if more than SF = 7 


is used then more than one card will be needed for each of the input 
variables. Here we will list the variable arrays as if only one card will 
be required for the input. 
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ALL FOLLOWING 

CARDS, FORMAT (7F10. 2) 


FCD 

array of frequency values 

CHZ) 

CORR(I) 

an arbitrary spectral shaping factor, 
normally set equal to 1 

(-) 

XLL(I) 

correlation length in the stream direc- 
tion, based upon the longitudinal 
component of turbulence 

(CM) 

YLL(I) 

correlation length in the direction across 
the wake, based upon the longitudinal 
component of turbulence 

(CM) 

ZLL(I) 

correlation length in the span direction, 
based upon the longitudinal component of 
turbulence 

(CM) 

XLT(I) 

correlation length in the stream direction, 
based upon the transverse component of 
turbulence 

(CM) 

YLT(I) 

correlation length in the direction across 
the wake, based upon the transverse com- 
ponent of turbulence 

(CM) 

ZLT(I) 

correlation length in the span direction, 
based upon the transverse component of 
turbulence 

(CM) 

The following 

variables, UlUl, UlU2, U2U1, and U2U2, are spectral levels 

given in dB relative to the RMS levels for the component. The 

result of 

this procedure will be: if these turbulence quantities are given in dU 

down from RMS 

on a 1/3 octave spectrum then the radiated sound 

will also 

be 1/3 octave, 

, If 6% spectra are used then the radiated sound will be 6Z 

also. 



UlUl(I) 

longitudinal component of turbulence 

(dB) 

U1U2(I) 

shear stress component 

(dB) 

U2U1(I) 

shear stress component 

(dB) 

U2U2(I) 

transverse component of turbulence 

(dB) 
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All input data are printed out in the physical units used in the 
computations . 


3.2 Subroutine WAKE 

This subroutine computes the wake related sound generation by the 
rotor. The model used to relate the force fluctuations on the rotor 
blading to aerodynamic variables is the single airfoil model. The 
expression which results and is here computed is 


p<x,(u) 


2 


1 

16 tt^ r^ 


I 

n=“«j 


B 


I 

k=l 


(^)2 

a 


^ -wk sin 4* 

d^(S^, £-nn) ce^, f-n« . 

o 

For the present case, the force fluctuations are 

A 

f-nfl) d^(6^, f-nfi) « S 6^ Cw-lrnfl)^ 

3 3 ^ 

X 

£=1 m“l 

where: S = span, 6^ = wake thickness, = correlation lengths (i «= 1,3), 

^l^ro " spectral level of turbulent stress, w = rad"? an frequency, U = rotor 
speed, and f = frequency in Hz. is a Bessel function of the first kind 
of integer order. 

The output from this subroutine is the radiated sound pressure level 
expressed in dB relative to a reference pressure of 20 pN/m^. 

3.3 Subroutine ORTIll 

This subroutine computes the components of the position vector of the 
sound observer in the coordinate system relative to each of the rotor blades. 


The product of these components with the force components on the blades 
yields the force in the direction of the observer. To do this the position 
vector must first be expressed in terms of the rotor coordinate system. 

It must then be rotated by the blade stagger angle to be in the blade 
coordinate system. 

The observer position vector components for a typical blade now may 
be written: 


R 

r. ~ (r sin - R cos 0)/[r(l ^ sin (}* cqs9) ] 

R 

r« = - ' sin 9/[r(l — — sin (j> cos 0)] 
b s 3r 

R 

g 

r^ = cos (|i/(l — — sin cos 0) 


The change from the observer coordinates to those of the rotor blade can 
be expressed in the rotation 


sin0 cos 0 


0 0 
\ cos 9 sin0 


n\ 

1 

0 


and the further rotation by the stagger angle, x. 




The stagger angle, X, is the angle between the airfoil chord and a radial 
plane passing through the axis of the rotor. 



3 , 4 Subroutine TOHE 


Subroutine TONE as currently set up computes the sound generation 
related to the steady state blade forces. This is the ‘'Gutin" mechanism. 
The program has been developed for a rotating steady state stress term also 
but as yet it is not known how to model this in terms of aerodynamic 
design parameters. 

The expression which is computed for the radiated sound pressure is 




(liR sin 6 
s ^ 


■) 


where: B = number of blades in the rotor 

F = blade force in the direction of thetbserver 
r 

R “ source of radius 
s 

(j) - observer angle to rotor axis 

r = distance to the observer 
w = radian frequency 


3i5 Subroutine 0RTH2 

0RT1I2 does exactly the same computation to the blade force as ORTHl 
did for the random blade forces. The procedures are parallel. 


3.6 Program. Listing 
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c 


L 




L 

C 


c 


c 

I 

t 

c 

c 

V. 

c 

c 

V. 

c 

c 

c 

c 

c 

c 

u 

c 

(_ 

c 

c 

c 

c 


h'KuGKAr''! i^uT jk( Ii^KUT »T ~ Ii\Pwl » IAPEI£j = uuIHuT ) 

rti.RUi-'YMrti'U C 6Gur\L/ UKi'ibkA i 1 uR I" KOV' F 1 dL AuL<-' i<l/TORo 
• •oiric. RAi^iATLlj iiQUfviu oPubrK*j)'i i b PRLtJ 1 bl LiJ 

OM bTc.rtu>Y LO'^t^S mNu AiRtOiL wnrsi;. olATj.il iv.rtL PROPc-H I i c-'J • * • 

Lili'ihl'iSIUN LilUl{5o) »uiU2(^'-') f\J'2.U7(b^) »F(b'J)»FT(2'JO} t 

iKhliT (200) »r<rlCoi3o) »CuKk(o-) >Ai_l(3^) ?aLT (ju) »YCl(3^) »VC.F(3o) 

^/^lL Tllu -JuuKLiuliMo lO iiNPul mV LuRuifi(j,\ I ht, 

Cm-i- uAlMiN(U4,ui»Ul'J<i-»'j^'-'i*'j‘^'>*i>F»KAo»i’(x»09(-rl»Ai_Ai'i»KO»F i''iA A > F tj » a1* > 
UL,yL»XLL»XUT»VCi-,YC FjCoaa) 

CAlL rtAAh.(F»>-'lui»Ulu2>sj2oijoi:o2»iMFsKPo>Ft'iAA9i’;o*tM'i(Jo»a» 

1 » XU Ai'i >Oi)ZL»Yi-»XLi-»Xi-T »VC>-*vfCF » ) 

,Vn<L v.OF'F-’'-' fl-O oRuAl. oUO/V oLjVRa I I ON 

UALu T vJi'i Lli\i»L;»ru9 Kriu 1 * AL/-M-. j rsf-'o » » u i » r f ) 

lOKii UOriPoTLS Lt'/'”Ub OF HOitu 1 ul\L ui_(si_Km 1 i,u ljY i\uiMTiNo 

oTt-ziUY ljU'^OL rOR'-i_ '■'Nu' i*.L/\N v/sU<->L-j <jr oIkUoj 1 lNoGR 

OF Jp 

FNO 

-^w«LiAooT INl. uA I M O'! { ui'-'l j oi(j 2 »^j 2^2 ’F ♦i\Pwi»f 4 l >q »i./ 1 *AuAi'. »i^o » 

IFi'iAX 9 F ri » I'lF s^LjYljXLl jaLI 9 VCi-»vCi » ) 

R[_hi>o ano oOiNUiiiON'^ OMiM 

DI^'iLi'- I'Ji'i UlUl(bw)s01Li2(bL,)»o2oi(ow)*u2'-'2(b^-'l»F(Fm) 

DI--.[£noIoN YC(30)»r0n;f\(3u)»XLL(30)>(’i-L(3-)»ZLL(3v>)9 
2 XLT(j«jj »V^I (3o) 9ALT(3uj *vC'I li«)9VUL(3-') 

Rp;-: = roToK opulij (kpm) 

Pi ” OO^Oi^VlR RmuIuo (iY) 

u - NuMljuK of RC'I^R dLAulo 

oi = oLmiJH RLLAFiYc: UtLOUllY l|Vi/oLv.) 

XL^Iv - O I i^OLitK A'-|GLU — F'ILmoOKOU FROF FhN AXIo ( i/LO ) 

Rb = bOURCt. Ka[) 1U3 t r. ) 

FA, AX = AX I .10, -I FAtrOoENCY oF liNTERtbT (HZ) 

Ft) = dLAoU FORCt. (X,b) 

aiuiti = R,.'lo Lo.NU 1 TUol NAl, TuKbuUUNUc - NurM-iAu 1 Ztu .vlfh 01 
uiotLi = Ri^io ontAi' — NOkH^LiodL/ >» I I M ul 
UtuiD = R.'ib bnEAK - NGKi'i'.'L I olu wild ol 

uZJ<ic = KHo FkAt'io\/EkoE Toi^tSuLLi'n-E ~ iM\jKi-,AL 1 XUU .* I 1 ti ul 
ZL = bP.AN UlhUNbloN ( C,"i ) 

ZL = WAXU TIICKNE.«;S (CM) 

CHD = BLAlJE CHORo (C.s) 

NF = NOfviDuR OF Frfq^<-NCY POiNFo 
F(I) = FRlwOENCY array (HZ) 

'-0RR( i ) = rt oPliuIrmL viiH'jriTii\u rUN^TlON iNOKPiAULY uRilr) 

XLL(I) = .jTkUAFi UCRkULAT luN Lti^oTri F Uk ululO 

rULdl - (MORHmL CoRRi-LATION Lc.NuTii for i^lUlu (UM) 

ZLL(1) = OPAN UORRELMTiUN LUNbTn l-Uk Uloio ( C,’^ ) 
iLUfCl) = oTRLAM UORRtiLATlUN LUNt-.Tn FOR o2U2o (U;,)) 

YUT(I) = normal CoRRuLATiOiN LtNuTd FOR U2U2 tUM) 

ZlT(I) = oPAN correlation LLNuiH rOR vj2u2 (CM) 

reproducibility of the 
ORIGINAL PAGE IS POOR 


n n o n 


21 


Ulul(I) = DlJ down r ROi'i L'lUlb I Ri'i^ ) 
= DB DOrtiN FH10|V( UlU2d 
u2ulll) = DB DO*'iM FROM O/iUlti 1 KM^ ) 
u2u2{I] = DB DOvviN FROM U2U2B (isMb) 

DATA AO,Co»Pi/3A2. 


1 


2 


READ { 3 s i ) KPM s Rl > R » iJi » XL At'i s Ko » Fi'iAX 

RLADI^jil Fosulolf^joiu d»j 2 ^jll)»o 2 ^ 2 d 

f«EADl6»i) ZLjYL»XL»CHD 

CMFT = .0328 

Ri = RHfCMFH^luO. 

ul=Uli^CMFT* 100 . 

KS = R 5 ^'C.MFT^fl 00 . 

FB=^Fri«-2.2i 

AL=ZL!^CMFI 

Y u=vls-c:mf 1 

CriD = CMD«C.’lh T 
format ( 7F 2 ) 

TPI= 2 .JtRl 

XLA.-'i=XLAi-.^RI/l 8 J. 


kP j = TPI4fRt^(v,/60. 


READ(5»2) NF 


FOK.mAT{1I3) 
k £ AD ( 5 » 1 ) 
.<EAD(b»l) 

I'i lAD ( 5 » 1 ) 
RtADl 5»1) 
K£AD( 5»n 
READ( 5? 1) 
Rt^AiKS 5 n 
READ( 5» 1) 
i<£a\D t b 1 1 ) 
k t AD ( 5 o 1 ] 
READ(5»1) 
RbAD( b » 1 ) 


IF« n » I = lsNF 1 
iCORS^l I ) » I = 1»RF ) 
(XlLC 1 ) » I=1?NF ) 

( VLLt I ) , 1=1 »kr 1 
!ZlL( I 1 » I =] »HF ) 
UlHI ) ) I = 1»NF1 
i YlT( n * I = 1»KF) 
(ZLT ( 1 1 > I=1»NF) 
(olollj)»I=l»kF) 
(U 1 U 2 ( 1 ) » l = l»iMr ) 
I o 201 1 I ) » 1 = 1 »tvF ) 
(02Li2( I 1 » l = l9NF 1 


C IiViTRUDUCE SCALIko FACTUKB 

oCL=(CHD/ 3 .Vi)*^^o, 8 *t 01 /lbv.. )’f^'-(-o, 2 ) 
oCF=( 01/ iPLf, )**1. 2«- 1 CHD/3.^ )** t “U.8) 
F.'1AX = FMAX*^CF 
XL=XL*bCL 
VL=YL*oCL 
ZL=ZL*bCL 
DO 30 I=1*NF 
F( n=F( I l*i)CF 
ALU I )=XLLt 1 )*SCL 
YlL( 1 )=YLLl I )*5CL 
ZLLl I ) =ZLL( 1 )*jCL 
XLT C I )=XU I U*SCL 
YLl {. 1 )=YLT ( I )^^SCL 
30 ZlT< I)=ZL1 I I)*SCL 
DO 20 I=1»NF 

VCL( I )=8»*XLLl I }*YLL( I )*ALLt I )/2B35^» 
20 VCT U )=8.*XLT( 1)*YLT ( n^ALFt l ; /2835^, 
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C rtKITE Tit INPUT UA-TA 

3 PORiviAT { ihi » 54t't AtROuVNAHiL ijUUNtJ uLNtKaTiON i>/,uLT 1 dLAdtU kOTOR 

10//23H OoohRVtR i^lSTaNUt. ( i-l) 1 F 7 . 2 » 1 > 1 2n RO 1 OR SptEO 

2 iri 0 «<i/i^n STtAUV tiLnuE hORLt; ^ F-i- i- • s i«J a » uLAUE dTAGOER ANOLt 
3 if i 0 * 2 / 17 ti . NUMbLR o>‘ ulaix-o ir 3 » v » 2 'jA 9 14 n oO'-'R'-c. RAuxud 

^-xF'J .4-/ (ia-H oLAUt RELmIIVE VcLoh I Y it*o» J. j iuA »lLin i^AX, fKtvxLrLNLY 
5il-iy*2/lwn Uiul(RMi)} iExw,d»L'A« i-rl Jiu^lRMo) lLiva,3»PA» 

6iuH IE 10 • 3 » 3X » i ori iblo»3////) 

WRITE(6»9) XL»YL 

9 FuR.'iATUlri v^AkE VoLU:»iP EtZ) = lE15»bj3Xt 
IbH L( Y) = IfclS.'i// ) / 

VYRITE(6»1H 

11 FORPiATtluOH frequency LtZ) LlY) 

1 LU) ji-aUR FALtOK 1 

rtRilLt6»l<!i) Iptl) »^LLl I}»YLLii )»Xll( I )»CUKl\( 1)»1 = 1 jNF) 

(vRITL{ 6»1^-) lF(l)»ALrtl)»YLlli)sXLl ll)»COKi\(i)»i = l>Nh) 

12 format (5E^O, 5) 

R I T E ( 6 » 4 ) 

4 FORcATdAH bPEClKiJM UlUi //) 
wRITElOsP) (F( I ) »U1^H 1) » 1 = 1.NF) 

f}F0RPAT(4(-LFl'j*3JlF2'-'»4)) 

•v K I r E C 6 » 6 ) 

6 F'-'r<*'tAl t / // 14H jRErTrSwn ijiuZ //) 

«RlTE(6»3) (F(l) »Ulu2tl)»l = l»t'*F) 

rtRlTEC6»71 

7 FUR.-iATI ///14H oPErTRuM u2ul //) 

.mRI 1L(6*3) tF( 1 } ,U2 uK I ) »l = i»NFl 
rtRlTEl6»8) 

8 format (///14H EPEcTRoM U2^2 //) 

*RlTEC69b) (ptl )»U2^2( 1) *1=1 »NF ) 

C CON'-'lliON I ORlJuLi- NCu. pa ) a to LxNLaR ooaLu 

DO 10 I=:1»NF 

XLLt I )=XLL( I 5 / t 2» A4*12. ) 

XLT ( I) =XL I ( I ) / {2* <^4^12- ) 
ulult 1 3 =uX 01 b/ luloH n / 2 <^* 1 ) 

ulu 2 ( 1 ) =uiu 2 b/ Clo,**iui«j-i(n/ 2 o»)) 
u 2 ul( l)=u<xulij/ti‘j»**( 02 'jlll)^ 2 'j*} 1 

10 o2u2( I) =Uc026/ t lu.#xM02u2( I ) /2^^. ) ) 

RETURN 

END 

T INt (^jAnEIF »Ulul»Uio2 >'j2'-'X»o2u2»NF»KPo9 Ki'iAa » Aij » rs.H'^o » B » 
lr(l jXLAJ-nui »ZL jYL jXLL »XlT »vCi-» YCl »Cui\K) 

C Cui'iPUTEb oUUNU RADIaTION l-Rui-i rta|(,L I uRouLEinOL PKcPERT i to 

U li'iti\5 Ion F t3^1*UlUl(3^) sulutldo) »u2ol(3'^3 » U2U2 I 3xi ) » B£b ( lu Ju ) » 
iBENlluuu )»PhH2^}*5X13v»31»j1'3o)» 

2RHOd(30l »V1(33»VC(3^) » CORK 1 3^) » 

3XLu( 3U ) *XlT ( 3U) *VcL ( 30 ) »VcT O'XJ 

UaTA CG9 AO*PURmO/1- 0,342. 14159*1.175/ 


23 


PRc,F-2U» E“06/ A0**2 

PREK = PREF^tPREF 

Cl = l. / t 4.«-p )*«-2 

TPI=2.*PI 

C2 = Cl^t RH0/CG)»*2 

REV^RPS/TPl 

AroT = Z.L*-YL 

c cui-ipuTE tme Radiation kattei^i'i 

DO lOOU K.|<.= 1,3 

PHnNA)=4P*-!^FL0AHKA-l) 

PHIP=PHI ( ^A)^PI/l80. 

c co/ipDTE r,ic spectpua 

DO li I=i»NF 

DVL = VCH I ATuT^tCOHKl i ) 

DVr=VCT ( I J^ATuTs-COKRt I ) 

DVL=oORT C DVL) 

D\/T=iQRl (DVT) 

DV=Q. 0 

oX( I >1 ) =-uioi ( I )*d''L*oOR I (XLLl I ) ) 

<=X(l»2)='-'2u2tI) *DVT*ouk r I AL M I 1 ) 

DO 40 J=i>3 

40 rtPlThtOjAl) ^X(IjJ) 

41 F ORr'i A TllUfl SXll»J) = 

DO 32 J=l>3 

32 VI {J ) =5X( I » J) 

V3=ulu2( I ) 

V4=U2U1 ( I ) 

FRA=F ( I ) *1 PI 

CAi-E OKTH-^(Vl»FiAA»V39V4»PHIt^>Ai_Mi't»£}»«l»i\o»DDl 

11 41(1 )=DD 

DO 2000 K=1,NF 
C3=(TPI*F( A) / A0)**2*-01**2 

C TrlE aRGUEi''iENT FUiA Tt-tti bEobEL hUNClliJN 

ARo = -TPI^^r ( a)*Kj*SI NtPHIP)/Ao 

C N -|mX 14 T I U- ^InXlt-'iUM NECtobARY \/riLot.b (jt" i') ruR Ii'u DLbot-L FoNO llON 

N-IAX=IFIX(Ff'iAX*TPl/RPb)^t2 + l 

in;-:ax=25o 

IFtNi'iAX.Gl.Auo) nNAX = 4^o 
i F ( ARo • NE • • u ) uo TO 6u 
bE4( i ) =1«U 
DO 61 LL=2»NKAX 
61 BEStLL)=u.U 
GO TO 62 
60 CONTINUE 

DO 63 LL=i*NMAX 
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63 ti£i> ( LL ) e tj 

CALL NBESJt AKG»NfiAX»bE:6) 

62 COMTIiMuE 

C ^jQuARE THc bEo^EL F'^lMCTIU.No 

00 44 J==1>NMAX 

B£N( )**( J-T )*bE5( J) 

BE-Jt J) 

44 bEM ( J : = bbi'lt J ] ( J ) 

C PROVIDE Tnt 6UMiv|A.M05 

bUMl=0,0 

DO 14 N=i»NMAX 

FP = MK)+FL0ATIN-1 )*KL\/ 

IF (rP.GT.F.vAX) 00 TO 13 
D‘ji‘*l = r3U'-<itFP»F>'Jl »lsNF) 

Doi'U = DUMlit(TPI*FP )**2 
GO To 14 

13 DOrtl=0,0 

14 6u,*il = SJPIl+D0[''ll^nES(N) 

DO 12 N = 2»N;'^AX 

FN] = F (K) -float lN-1 1*REV 
IF(Fi\.LT.-Fl'iAX) GO lO 13 
Fi\ = ABO(FN) 

1= 1 tJLu 1 1 FN »F’0] »1» NF ) 

Du,^l = DU:vU*-( TPI*FN1 4t-^2 
GO TO 12 

15 DUM1=0.0 

12 ooMl = bUM14-DuMl*BEMlK) 
wRITE( 6»31) bui'll 

31 FORiMATtbH ooMl lE2^.6»3n D'-'Mi 1E2^*5) 

C THIo COMPLETc:) THE DlPOLt 1 LKi'i 

c V»RITE The i^iODoLu5 of the KADlAlcD oUUWD DEHolTY 

VARi=C2*CJ*SUMl/PpEF 
rtRlTE(6»50) F<X),VAH1 
50 FOR.'iAT ( lFiL,3»lE20,5 I 
2000 RH06 ( K, ) =C2*C3*S0 m 1 
VaRITE{6»l) PHIIKK) 

1 FOKi*.A i ( iH-L » 37H oPecThAL DI o l K I doT Il-N uF r^uTuK ooJHD /// 

130H mrAKL TURddLLMCE RELmIilD iOOND //17n RAui/vliON ANOli- ir7.2// 
210X»i0H FREQUENCY lOX.tiri bPLlDb) /) 

DO lO I=lsNF 
VARl=RHObl I ) /PREF 
IF! VARI.lT.1,0) go To lu 
KH0b{ U=1U,»ALOG10( VARl ) 

10 CONTINUE 

WRITU692) (F(I) ,RHU6( I ) »I=1»NF) 

2 FORi-iAT (12Xf lF10.2,lUX»iE2u.4) 

1000 CONTINUE 

RETURN 

END 
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ouBkOuT Irtc OKlHitSX»Kt-*^»ulu^»Uiiwl jPHI » XL Ai''i »B »k1 » »0KD« ) 

L '-OhPuTc.6 Trill ooM of" T ML LOnPOt^ii-N fi3 IN ThL iJiRuLliON Ol" i '^L 

L OLiiiRVciR rRoM LALH liLAlJU Of Tnc. ROTOR 

L)I.»iENi.ION ZY{3»3)»TmL{3^)jYX(3»J)*oX13)»^TY(3»3)»i-4(31*V(3) 

C iNTROL>uLn TmL COMpONi-NTi OF 1 Fil^ST ROTATiUN 

PI=3.1A15V 

ZYi 1»1) =blNCXLAM) 

ZY( l,2)=COS(XLAM) 

ZY{ 1 j3)=0.U 

ZY(2>n=-CUS{XLAN0 

ZY(2»25=5IN(XLAM1 

ZY(2i3)=0«L 

ZY ( 3»l5=ue'J 

ZY ( 3 »2 3 =U#U 

ZYt3»3)==l*u 

^ TiIl NLXT rotation Oi-PUROb on 1 ni_ iR^TinL ljlAi^L PObiTlON Ab bOi-o 

C THc ObbCRVtR RAuiiUS VhCTOR 

N D= I F I X ( 3 ) 

DRDR=u# j 

00 10 K = lfNtJ 

THh( O = 360 ./d*FLUAT (k) 

ThFT = THE tN)*Pl/iaOo 

0:iN = Rl»l 1 •-Rb/pl^^5I^! I PHi )^tCUo ( Trtt'T ) } 

Rtl) = (Rl*-3lN(PHll -Rb*roiS t THET ) ) / 0£N 

R( 2 ) = (-Rb*oINnilcT) 1/DFN 

R ( 3 ) = ( RlftCOSlPHl ) ) /DEN 

ifX{ 1»1) =s~^lNl ThET ) 

rX{ 1»2 )=CvJo{THET ) 

irx( i53i = L/«'-' 

VX ( 2ff 11 =0.0 

YX t 2 »21 =U.u 
YX(2, 31=1.0 
YX(3>1 1=C03{THET1 
YX( 3»21=01N(THET) 

YX( 3»31=0.o 

L NOW bUM Ol\ The. CON'PONLNTo To uqT rj.i\aL VaLolo IN ObuLKYi-R l-Rrt,v|u 

DD=0.0 
DO 11 I=l»3 
V( 11=0.0 
DO 11 J=l*3 
DO 11 L=l»3 

11 V{ I )=R{ LI^ZYI 1 . JMS-Ya IJ?L)+Yt I ) 
wKITE(6»3b) {VlIl*I"l*3l 
30 FORi-IATtSH V{Bl 3^15.3) 

C buTLRMlNt Tht COMPONc-NTb IN T ML DlRtCfiON 01“ THE OObtRVtR 

C DIPOLE TERMS 


n n 
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00 I=l»2 
DO 12 J=l»2 

IF( IotQ*l*AND.JoL0,2 ) uO TO 13 

1 F ( J • t.Q* 1 • A(S)D e I » t Qo 2 ) OO To 13 i 
VA'R 1 = SX( n^^SX( J) 

GO TO 12 
13 VaR 1 = -U1U2*^'^2 
GO ro 12 

131 VARl = -U2oi'*Ht2 
12 DO = DD+V{ i )*VAHl*v ( J) 

10 DRDR: 'RDR+OD 

wRITE(6»20» DRDR 
20 FORMAT ( 1E2U. 5 ) 

C Tnii. C0 MPLlT:i: 3 CO^POTATiON Qt" TnL i-OMPOMc-N f ^ Tnc 

C uiRi_CTiOM OT THc. QbOtiR\/LR 

RETURN 
END 

bJBKuUT INE ToNEt'^l » B » FB » i^HO » X LA.-i » nPb » Kb » ui »F ) 

COMPOTES Pure tones OENEKATlU oV bTcAuT uLAOE Fv,'KCES ANo 
Thu MlAN vAL>bES OF ThE biRtOb I UN-jCJK 

0 1 riLi'^ b I uh btb(2“^ J >bX(3) 13»3J » » lYY >F( 2'-J'^) > 

1RHO(2UO) »HHI{20)»TT(3u.) tD0<3^! 

OAlA AO I P ■» » CG/ 3A^i o » 3 • 1 A i b y i 1 ■ u/ 

PREF=20«E-06/A0**? 

PREF=PREF"PREF 
RATM=1,173 
RATM=0*0 
NI3=IFIX(B} 

0 PROViou fiiu TuRoOR ANO lOKOL ^-OhpONoNTb 

C1 = 1./(2««PI*A0**2^^R1 )-«-*2 
RLV = RPS/ t A,^fPI ) 

FAX=:0.05*FB 
FT = 0,95«-Fb 
DO I 1=173 

1 SX( I )=0,U 
SX( 1 )=-FAX 
bXt2 )=FT 
DO 2 1=173 
DO 2 J=l73 

2 oYYt I 7j)=0.u 
SYYdf 1)=1«U 
DO lOOo K.K = 1?1 

PHI (KK)=30. + lu.*FL0AT{tCK.) 

PHI ( tUiO =6U. + 1 U**FlOAT tKK) 

PHI ( KK) =80. 

PHIP = PHI (««K)*PI/1rO. 

C2=(3/RS)**2 

C3 = tB*RATM^tUl»Ul/ {A0*RS*SIN(PHIP ) ) )«-*2 


nr' r. n 
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C compute TriE OPEC r RUM 

caul OKTH'^lOXs6YV»PHlP»ALAi''l*i}>''4i>rv,b»DD»I 1 ) 

C COMPUTE 5 harmonics OE dLAOl prtoSH^JL hRLQUtNCY 

NH=10 

DO iO N=1»NH 
F{N)=FLOAllN)«REv‘^tB 
ARG = 2.i^PI*E (N)*R5^{-SirM(PhlP) / (^U 
NN=IFIX(B)*N+1 
call NBESJ ( ARGjNi'i ,BES) 
wf?lTE{6»30) ARG»BESiiMN) 

30 FURMATtl6n i3ESSEL FUNCTION 2E2U.5) 

TYY(N)=0.0 
TX { N ) =U o L 
DU AO X=1»NB 
DO AO L = l*NcJ 

TXt N ) = 1 XI N ) +BES (NM ] *LiES t i’lN 1 ^ ) ^^dDI L ) 

AO T YY ( N ) = TYY t M )+BEo (NN ) ^dE’^J I NN ) ’^T i ( k )*1 I I L ) ^^2 • *P I *F ( N ) 

Tnlj COMP Lc Til Tf'E oPECTRrtU LlVlL '-UMPu Tm T i UN - NOrt cAlL ROuTiNL 
T o PLRFORi- COORlJi MATii- iRr^Nj^ORM^T iONS 

VAKi=Cl*CA*TX(N)/pREF 
VAR2 = Cl*CJ»^tTYY(N) /PKLF 
.'vrti TE I 6»2U ) FIH) , VAK1» VAK2 
20 FOR.'iA T I 3 ( IF I'J.Z >2rlb .3 ) ) 

FiNN = FlOAT ( N ] ^thLCAT IN ) 

10 RHU( N ) =C1*C2*FNN’'‘'TX ( N ) 

ARirE(6»5) 

3 FOHI-iAT ( IH J. j 37ri ^PfCTRAL i-' I T K i LuT i UN uF koTuK ouuND ///) 

DO ii I=i»NH 
VAR<i = RHO( 1 ) /PKEF 
I F I VAR2. L I o 1 • vj ) pO To 11 
RH0< I )=1u,*ALOG1^ ( YAR2 ) 

11 CONTINUE 
WRITE I 6 *3) PHItr^K) 

3 FORriAT ( / // / /36H ^PEC ( RAL D 1 f Rl oU T 1 UN TciNE o'jJhD /// 

117H RrtDIAliON ANclE 1F3»2/// 

29H FKEQ(HA) 11X,8H SPUDcJ) //} 

wRllE(6*A) (F(IJ)RhU(I)»l = i»!^lril 
A FORriAT ( A ( IFI'^ *2 » 1P2>J • 3 ) 1 
iOuO CONTINUE 
RETURN 
END 

::iUBKOUT INE OK Th 2 I 5X * b Y Y » PH I »XLAH,a 5 Kl »Kb »DKDR » T «TR ) 

CQMPUTEb Tht bUM oF TtiL '-0>MPUNi-NTb IN T.iL uiRi-CllON OF 1 rlL 
OUERVlR rRO.y EmlH oLAl>l of Thl rotor 

Dli'iLl'iblUN AY(393) 9 THl(3o)»Ya 13»3) »oXl3) ♦■SYY(3»3 )»k( 3) »V(3) 
DIMENblON DRDRI3U1 ,TRTR(3i^) 

c Introduce thc componunts of the first rotation 
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pi=a. 1415V 

^Y ( 1 »1 ) =oiN(XLAM) 

ZY ( 1 »2 ) =C05(XLAi''i) 

XYfl»3)=0«U 

XY<2»1)=-'C0S(XLAH) 

ZY(2 ?2i =SIN(XLAM) 

XY{253)=0.U 
XY(3»1)=0,U 
^lY{3»2)=u.u 
2Y( J»3 ) =l*u 

C Tut NtXT t^OTAnOlN Utf-ENDi) ON Tnt iNiUAL dLAuh HOSIiiTlON A6 dOtb 

C the observer radius VECT'^R 

NB=IFIX(B) 

DO 10 K=1»NB 

TH£( K)=36u./t3*FLOAT<K) 

rHET=THE ( \)*PI/ldO, 

DE,v)=Kl*( 16-KS/R1*SI<^ tPHi )*C^ol InEl ) ) 

RI l)=(Rl^oIN(PHl ) -RS*CUd ( THE f ) )/DtN 
R(2)=(“R^*SIN(THET) ) /DEN 
R( 3) = (R1^^C0S(PHI ) }/DEN 
YXn ,1 )=-SlN( THET ) 

YX( 1 »2}=COi( THET) 

YX( 1 »3) 

YX {<i s 1 ) =iO«U 
YX(292)=U«P 
YX( 2 * 3 ) = 1.0 
YX<39l)=Cod(THET) 

YX(3»2)=SIN<THET) 

YX ( 3 *3) =0.0 

jNOo Su‘'1 on THl EO‘/PGNtNlo lu oc I riN^L VmLuco iN Qu^cRVl-R 

DD=0.0 
TT=0.0 
DU 11 1=193 

vn ) =0.0 

DO 11 J=ls3 
DO 11 L=l93 

11 V( I)-k(L J*YX( J»L)*ZY ( I »J )+V( i ) 
rtRlTE(6i.2U} {V(1)*I = 1*3) 

20 FuRi'.mT [ //31ti ObSERVtR VtuTUrt IN dUMUE bAdlb /3t-2o.5//) 

C UtTtRMlNt^ Tnt COMpONtNTS IN Tfic DlRtCTiOM OF Tni_ OoiiERVtR 

C DIPOLE TERFiS 

DO 12 I=l>2 

12 DD=D0+V( I }*SX{ I ) 

WRITE(693U) OD 

30 FORMAT tl.Sn FORCE COMPONENT iE2u.5j 

C OUADRAPOLtS 

DO 18 1=1*2 
DO 18 J=l»2 

18 TT=TT+V( I )*SYY( I.J)*V( J) 

DRDR( K)=DO 
TO TRTR(K)=TT 

C Tnib COMPLc.TtS COMPUTATION OF iNt COMPONENTS IN THE 

C UlRtCTlON OF THt obi^tf VtK 

return 

END 
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4.0 PROGRAM TONE 

This model was set up in a computer program so that direct numerical 
computations of the time series described by the equation for the radiated 
sound could be made. This direct approach of calculating the actual time 
series is the difference between the present approach and those of other 
rotor models. Typically the analytical results for the rotor are Fourier 
transformed and reduced to a closed form (such as program ROTOR described 
herein) . This approach normally requires that simplifications be made so 
that the closed form can be obtaired. The advantage of the present 
approach is no simplifications In the analysis need be made. 

The time series is computed for a sufficient number of values so that 
a good statistical sample will result. The time series is then Fourier 
transformed using a Fast Fourier Transform algorithm to yield the spectrr,.' 
of the radiated sound. In other words, the model is used to generate a 
time series which is handled exactly as one would handle experimental data. 
Very complex source histories can be studied with no more effort than the 
very simple steady force case. For example, rotor wobble can be simulated 
in order to study how it affects the generation of sound; the individual 
blades can be allotted to vibrate as might occur in flutter cases; or each 
blade can be assigned a different value of force to determine the effect 
of such blade to blade differences as might be caused by manufacturing 
imperfections. 

Of most importance, is that complex random effects can be included in 
this model, also. Such mechanisms as force fluctuations related to wake 
turbulence, inflow turbulence, or even rotating stalls can be studied. The 
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way this must be done is to express the force fluctuation in a tltne 
series that has the required statistical nature. If the spectrum of the 
force fluctuation relative to the blade is known, then it is possible to 
generate a random time series which wJ.ll give the same statistics as the 
knoTO result. The methods for generating an auto-regres .on analysis can be 
used to simulate the spectrum of the force on the blades. Once the time 
series are known, the results from the different blades may be combined In 
different ways to simulate different mechanisms of sound generation# For 
example, the effect of large scale turbulence can be simulated by allowing 
several of the blades to experience the same force fluctuation, or allow all 
of them to experience it. The case where there is no blade to blade corre- 
lation would model the wake related sound generation. 

The idea of the model is to provide a tool where the various mechanisms 
which might produce appreciable sound generation can be evaluated. In addi- 
tion, it is sought to set up the model in such a way that its usage will 
be reasonably simple and straightfort'/ard . As an example, if one wants to 
study the effect of different forces on each of the blades, he merely reads 
in the array of steady loads he wishes to rotor to have, and that is all. 
Rotor wobble is accounted for by replacing the velocity of the source, 
the rotor velocity, with a time dependent function describing the wobble 
characteristics. 

4.1 Tone Noise Generation Model 

The generation of pure tone sound due to the steady state forces on the 
rotor blading is studied in the model. In addition, an auto-regression 
method is used to write a time series which approximates a random force dis- 
tribution resulting in broad band noise generation. 
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The model Is based upon M.V. Lawson development for rotating point 


forces . The development is done for a force per unit volume. In this 
representation only the point force case will be considered. Extension to 
the distributed case can be treated as an array of point forces where the 
array is distributed along the span as well as circumferentially. This has 
been done by the author and the mechanics are straightforward but the 
expense in computer time is correspondingly increased. 

For a point force, the radiated sound In the far field can be written; 


x.«y. 3F. F 3M 

pC#.t) - t 

4iiar2Cl-H)2 ‘r 

o r 


where the evaluation of the source terms of the retarded time 


= t 


= time of sound generation 

= time of observation, and j'x - *y [ /a^ = propagation time. 
The position vector r = x-y can be written in the form 


r-^ (x^-y^) in rectangular cartesian coordinates. 

The Mach number is the component of the source Mach number in the direc- 
tion of the observer. It may be written 


M = M • r/ r| 
r s ' 


or in components as 


’’The Sound Field for Singularities in Motion," Proc. Roy. Soc., A, ^3^, 

559, 1965. 


S' 


1 
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«r" ^ H^(x,-y^)/|t|. 

i=l 

The components of the blade force F can be chosen to represent the 
thrust and torque loads on the individual blades. Using the blade stagger 
angle, 5, write 

= -F cos e 
= -F sin 5 sin 0 
F^ = F sin C cos 9* 

Take the time derivatives of these force components; 

3F^/3t = 0 

SF^/St = -Ffl sin C cos 0 
3F^/3t = -Ffl sin ^ sin 0. 

The components of the position of the observer are written in polar form as 

- r cos 1}) 

X2 = r sin ^ 



The location of the source is also written in polar form; 



^2 “ is ° 

yj = sli. e. 

Finally, the components of the source Mach number are 
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=: 0 


1 



“2 

='M 

sin 

H, 

= M 

cos 


The value of the source Mach number is 

M = R SJ/a . 
s s o 

The Mach number component may now be computed: 


3 

i=l 

Substitution yields 

M = (r^ cos (|))(0)/r + (r^ sin ^ - R cos 0) (-M sin 0)/r 
r X X s s 

o 

- -(r^/r) M sin ^ sin 0 + (R /rV^iT^os 0 sin 0, 

Is s s 

so that to the first order 

M = -M (r^/r) sin 4 sin 0. 
r si 

The time derivative of the Mach number becomes 
3M /3t = -M (rj/r) fi sin ^ cos 0. 

All the terms required for the computation of the pressure radiated are 
now complete. The radiated sound is 
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pCx,t) = [• 


4;ra r^(l-M 
o r 


sin c(j - R COS0 

{ ^ (“fl sin C cose) 


- R sin 0 
s 


(- SJsin K sin 6) + 


- 3M r_ cos ^ 
1 r r 1 


[■ 


1-M at “ r 
r 


(- cos 5) 


r. sin (|» - R cos 6 
+ ^ (- sin ^ sin 6) + — ^ 


- R sin 0 


sin ^ cos 0] }] . 


And, expanding: 

pC^,t) 


= [ 


p -r^ H sin 4* sin 5 cos 0 


4:ra r(l~M )^ 
o' r 


R n sin ? 
s 


3M -r- cos 4 cos ^ 

Ccos^e + Sln^e) + 3^ t ; 

r 


-r^ sin C sin t sin 9 


]}]. 


The solution reduces to the form 


pC^,t) = 


4iTa r^(l-M 
o r 


[n sin 5 (R^ - sin <j> cos 0) 


. 3M 

- ' y — ' Cr- cos (j) cos ^ + r, sin C sin 4i sin 0)]. 

1-M dt 1 i 

r 

This will give the radiated sound from each of the blades. The above 
expression is to be evaluated at the retarded time. This will require some 
further manipulations • 




where r = [x-yl 


square the distance r 


m Cx,-y.)^ " (r_ cos + (r- sin (p ~ R cos 8 )^ + (R sin 8 )^ 
X i i X s s 


= r? - 2 r- R sin <ji cos 8 + . 

11 s s 

Now, r_ » R , so it will be sufficient for the far field solution to 

^ -S 

write 

r^ = r? - 2 r- R sin 4 cos 8 
1 Is 

or, equivalently, 

= (r - R sin (ji cos 6 )^ - r 2 ^ j,pg2 0 ^ 

A 3 S 

For the far field solution we can therefore use 

r « r. “ R sin 4 cos 6 . 

Is ^ 

The retarded time can be written as 

= t. - + Rg/a^ sin 4 cos [9] , 

where the [ 8 ] is the retarded value of the blade location. 


This source position may be written as 


[ 8 ] = ntj^ + 0 ^ = fiCt - + - — sin 4 cos [ 0 ] + 8 ^ . 

Now in the solution for the radiated sound substitute t^^ for t and [ 8 ] for 6 
the radiated sound pressure ' -comes 


p„(x,t) = ;; [n sin C(R -r. sin 4 cos [ 6 ], ) 

^ 4Tra^r2(l-[M^])2 si k 

3M 

“ ( 1 - T m 3 ) ' C cos 4 + sin 4 sin 4 sin [6]^)} 
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where : 


with: 


r = sin iji cos [0]^ 

sin ^ sin [03j^ 

[aM^/8t]^ = -Mg J2 sin cos l0]^ s 

[0];^ = net - r^/a^) = cn R^/a_^) slu ♦ cas [6]^ + 

®k t =0 ~ location of kth source at = 0 . 


The total sound from all the rotor blades becomes 

B 

p(x,t) - I p, (x,t) . 

K=1 ^ 

4.2 How to Use Program TOME 

As was seen in the previous section, the model used to predict the 
sound radiation generates a discrete time history using the solution for a 
rotating array of forces, then itreats that history in the same way one 
would treat experimental data. 

We previously mentioned the flexibility of this approach. This flexi- 
bility is achieved in the program by changing the program listing, A 
multiple option approach could have been taken but this is wasteful of 
computer time and assumes one knows all the options which will be needed. 
The reason for this approach was to provide a flexible tool; however. It 
does depend upon the operator’s ability to change the Fortran statements 
used to generate the time series. This is easily done since the equations 
are short and the program is straightforward. The part of the program 
actually used to generate the time series is about 20 cards long. 
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The form of the program set up here is for the case where the force 
on each of the blades ds the following function of time 

F(t) = Fg (1 + C^x(t)) 

where F^ is the steady state blade force, is a dynamic loading coefficient, 
and x(t) is a time series simulating the unsteady force fluctuations on the 
blade. Depending upon the particular mechanism, x(t) can he chosen to 
simulate either wake related force fluctuations or the fluctuations caused 
by inflow turbulence interacting with the airfoil. Processes of this type 
Can be represented using autoregressive processes. This will be discussed 
in section 4,5 below. 

4.3 Program Inputs 

Input to the program is done using cards. Four cards are required to 
run the program in its present form, 

CARD 1, F0HHAT(1I5) 

NOPT = 0, no plots 

“ 1, plot the time history 

=2, plot the sound spectrum 

= 3, plot both the history and the spectrum 

ROPT controls the type of output desired. Two CALCOMP subroutines are 
contained in the program so that the time history for the process can be 
examined and/or the spectrum of the radiated sound can be plotted. If 
neither plot is desired then listings of the time history and the spectrum 
will be the only output. 


CARD 2, R0ISMATC7F10.2) 


BE 

bandwidth desired in the spectrum 

(Hs) 

FMAX 

maximum frequency in the spectrum 

(Hz) 

See Section 

4.4 for the implications of these selections. 


CARD 3, FOR1IAT(7F10.2) 


ZR 

number of rotor blades 

C--3 

XI 

blade stagger angle, ueasred from blade 
chord to radial plane passing through 
rotor axis 

(DEG) 

Fin 

radiation angle, measured from inlet 
rotor axis 

(DEG) 

RPM 

rotor speed 

(RPM) 

RS 

radial distance to source 

(M) 

R1 

radial distance to observer 

(M) 

P 

magnitude of force on each blade 

(KG) 

CARD 4, FORMAT C7F10. 2) 


CDYM 

dynamic force coefficient 

(— ) 

A1 

autoregression coefficient 

(-) 

A2 

autoregression coefficient 

C-) 


See Section 4.5 for a discussion of the selection of the autoregression 
coefficients. 

4.4 Tina Sampling 

The setting of FMAX and BE determines the sampling rate and the total 


number of required samples. Letting 
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we can illustrate how this isdsne. First, estimate the number of samples 
required as 

51 = 2 f /B . 
max e 

Since we will use a Fast Fourier Transform to find the spectrum, the number 
of samples must be of the form N* = 2^ where m is an integer. We choose m 
such that 

2^"^ < N £ 2^. 

Then compute = 2™ samples of the time series. The new bandwidth 
becomes 



f 

max 


/N' . 


The sampling rate for the time series is 


H = l/(2f ) , 

max 

4.5 Use of Autoregression Processes 

A second order autoregression process is used in the current form of 
the program to represent the force fluctuations. This process may be 
written 


x(t) = xCt-1) + H 2 x(t-2) + zCt) 

where and a^ are the autoregression coefficients (A1 and A2 in the pro- 
gram. Z(t) is a random number selected using a randum number generator. 

A time series is generated separately for each blade. If there is to 
be no correlation between different blades (at? for wake related forces) then 
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Z will be taken from a different random series for each blade. If correla- 
tion is desired between blades then the time series can be selected to 
reflect the correlation. Correlation would be expected to exist between 
blades for the mechanism of inflow turbulence. 

The main idea is to fit the time series to represent what is desired, 
then let the program compute the jresulting sound radiation. 

The spectrum for the second order process above can be solved for 
directly. The spectrum is 

G (f) = Ho^/Cl + a? + - 2a, (l-a„) cos 3rfH - 2a„ cos 4irfH) 

XX z 1 2 1 2 2 

- 1/2H 1 f 1 1/2H. 

Selection of a^ and ^2 gives the spectral shape to be expected from x(t) . 

A plot of a^ and a^ for dlffereiit spectral shapes is shown below to 
give an idea of the values of the coefficients needed for particular 
shapes . 



LOW FREOUENfY SPECTRUM 


c' 
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A peaked spectrum is what we would expect for a wake related sound genera- 
tion mechanism. A lot? frequency spectrum xrould be expected for the cases 
of inflow turbulence. The frequency where the peak will occur is givey by 

cos 2-rrf^H == (1 - . 

As may be expected, some experimentation will be required to fit the process 
to an actual experimental spectrum. 

To obtain more complicated power spectral shapes higher order processes 
may be needed. These can be handled as easily as the above on the computer, 
but the preliminary efforts to get a particular spectrum increase. 

Tor a good general discussion of the above method the reader Is referred 
to SPECTRAL AtlALYSIS AND ITS APPLICATIONS by G.M. Jenkins and D.G. Watts 
( 1969 ) . 

If more simply definable fluctuations are needed these can be used for 

i 

x(t). Tor example, inlet distortion or the effect of another blade row can 
be simulated by a Fourier series representation of the force fluctuation 
resulting from the distortion. A simple cosine with period equal to the 
stator spacing will give a good qualitative model of blade row Interaction 
noise. 

The model can also be extended to account for span distributions of the 
blade force. This is straightforward but the computer time required 
does Increase substantially. 


nnno n on onnrnrinnnnrnr n n 


PKJGKAf'i TuNEt I NPUT wUTHU'f » 1 At^E5= I NPUT » T AH E6 =oU T t^Ul J 


CUHPUTh buUND RAl) I A T I OiN FkOi'! a KOTA Ti, MG ARRAY 
OF POINT t-OKCtO 

COPtPLbX C 

OltnLlMGIUN PKADt 1100) jTHtPt 2u) *TMt (2o ) »TftFKl2o ) » IritRPI 20 ) 
OIrtLMGION C( 1U2A) »B t ^ » l^-i4 ) » 2:^6 J »H( 3) 25b 1 » AMP ( lu24) > 

IOPl. ( 1024) 9FREw( 512 ) ’ HkA( j. >^ 2 A ) » A ( lo24 ) »X { h s 1.^24 ) »pA TO ( 2o48 ) 
E0UIVALENC£(C9B ) 

KiiaR = naRMOMKO TO ot ANaLVAtiu 
NOPT = OP I ION 

lJl, = t-RcQuLNC-Y 1 uTri JLoIRiilJ 

IR = MoMdCR OF RuTOR BLAuLo 

XI = bTAOotR /iNoLE» MtabORLi-' i' ROM FmM AX 1 b (Ucb) 

Prii = RAO i AT I ON A^GLE - i‘i t AiURct-' fRUK FAN AX lb 
KPri = RGTo-K GPEEg CRPM) 

Ro = bO'JRCE RADIOS (h) 

R1 = OBoEKVER RAujU^ liM) 
r = l-OR-wti ON LMcri ROTOR uLauc ttsu) 

COYN = OYNmMK LOAOiMG - Pck bcNT JF TOTAL hOkC.- 
Ai = AOTO-REGREooION CCM^TAMi 
A2 = AUTO-REGRESbION CONSTANT 

XCJjI) = Al^XlJjl-l) + A2i^X(j9l-2) -4- Xtl)' 

4(1) = RAMUCiM variable 

DAlA Pt jAg/3c141i'99 34.;./ 


INP'-'f Ml •iUCn TriAT C HAb LENGTh M = 2**Ml 

READtbsi) NHAR»M0Pf 

1 FORMA7(2Ib) 

READ ( b 9 2 ) BE 

OPTION = 0 NO PLOTb 

upriuN = i PLOT time HibTuRt 

OPTION = 4 plot SPECTRUi'I 

OPTION =.3 plot both TIf-it nl^TORY ANU bPtiCIROM 

kEAD( b»2) ARsXI tPHi * HPi-; , Kb , k 1 » F 
kEAD(b»<;) CDYN»A1»A2 

2 FURMAT( 7 Fi^Jo 2 ) 

C CONUITION INPUT DATA 

CMFT = 3b 48 
Rb = Rb*C«^FT 
kl=ki*CMFT 
F=F#2o21 

FmAX=FLUAI (N riAR)*RPM/6B»*4K 
FMAX=lOOOOo 
N1=IFIX(2 .*FMAX/Bf) 

DO 75 I=i»lU 


repboducibiuty op top 
original Page is pooh 


43 




I 

IF(N2.GT.iMl) GO TO 71 jj 

7b col'll I )\uE i 

7i NoAiMP-Nid I 

1 = I j 

IHE MAXI^1^J^; ^U^;3ER OF POlMTo AUE (lu2A) I 

1 

bt = 2 b-h FMAa/FlOA I I MbANiP ) i 

Xoi'iFG = KP.'.5^2B^Pl/60. ji 

XIP=XI*PI/18j« ^ 

PnlP=PHl«PI /Ibj. I 

-1i'<PHI=OIi^(PHIP) I 

C!^bPnI =COo ( PHIP J i 

bIi'JXI = OlN(XIP) I 

COPXI=COotXIP ) I 

X.‘io=i^bii-XOi’'iEG/AO ^ 

;mzr=zp I 

r(ji\;_L UiN IhL KC T aK uL'^^i-o Lr i_ L ^ i_* Y H ^ i\ Fs L. L f ^ i-| 

o t L L. K 0 O i-< iv E R p R U C E G 0 

IX^Z»*2A+b I 

NRANO = NSA.'.P |j 

D.J 31 J = i »XZi< 

lX=IX+IU 

GALL RANO t I u » i X 9 1 » R R AKD »<- ) :;j 

X ( 0 > ) — 1 c l) :-i 

.'5i 

A{JfL)=loW 

00 I=3»RjA,-:p i 

31 X( J»I )=Ai*XU»I-l i+A2^'X( I ) 

WRITE INPuI data 

S 

<1 R I t L ( 6 9 2 u ) A R » A J » P(^i » F 9 K J » r M I » i< i 

20 COR,’-R'i (lnj-»biii RADlAfiOiX (JF jo'J.'ml.' ^ KUiMli.-^o mRR^Y f G P' l. o T 

MvJi'hUlR Of POIgR i>L^'i-'<- -jlfiuvJLp (^fyjoLi- 

2lFb.^/ldrt ROlOR bpLLiy (RPM) if 1 • 2 / 1 2< s qL'>l»i- rORGL irlv^.2/ I 

bJiri RADIAL LOLATiOM OF Pulivif FoKCE 1E;j«2/ s 

4I6H RA(Jl/\liOR ARull ll-b*2/iOi> Qb^uRVcR RAjiuo iFb*I///) 
wRITh(6»2^) FPAX»HE | 

22 F JRi'-iA T ( / / X tiri I'lAX 1 FKlOulimCY G 

1 Ir 3 » ibi I uARuiviDT” ( ri2 ) iCiOot.//) 

..RITL(6»2j) CDYR9Ai*A2 I 

2 3 F OK i-i/-' I ( 2 7n HAROori FGaCu F lov» [ un 1 I utMo // A 

137H XU) = Al*Xtf~l) + A2*X(|-2) + U ) / | 

2di:)n OYN^'M^G COLFrlCiLMl - xrxo.J / | 

3oHAi = J-Flo*39iOX»^ilA£i = xFiu*j /) f 

iRlilALlZt- ljLADl POoiliONo I 

Ou x0l=l»RZR 0 

01 = l I; 

10 the t I )=2,*PI/2R*tUI-l» J i 

CoMPUTti time GLRiES uAOLiv ON Rt-lAlRlNu c>uFFiLitiNT | 

information for 10 HARMONlOO I 


r C -. 
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C jOt SAMPLL at 20 sAHPLtS PtK uL^^UL i,PACIIMo 

C oL.r 1 i.ML i i\Li^Li\lI'JT 

FLi = F 

Dr = 0.b/FiMAX 

D.,' 200 K = l»,\iSAMP 

JK-< 

T=DH^(U<-i. ) 

L >-0;«ipijlL I M f X CiM r- KOPi I 'l>- ^>KKaY 

= U*U 

Du 300 J=liNZK 
DJ = J 

IF(\.ML.n GO TO ^Cjl 
iHtPt j) = Tii;:( j) 
vj J r u 3 J2 

3^1 TMlP( J)x=ini:KP( J) 

3 02 Ool'tTINLlG 

rHC>'ItJ) = Xui‘iFG*tT“Kl/Avj + i\o/A>j’^oIuPnI'*C'^‘jt TUKr^tu) ) j + ij-iFtU) 
ThtKPt J) = iHEPt J) 

A, ,K = -xv,o<voIKP(IHoI;N( i,-a J1 ) 

uX) iK = -X /.o«XOMLG^' ^ I ^jPfl I ’’ Cuu ( I MuK < J ) ) 

P=Rl-:'?oJt^lf^PHHtCOG I ifiFM J} ) 

F = F->J^ [ 1 , + OnYfM^^XlJ».< ) ) 

vA:< 1 = F / ( I ^ sK - « 2 ) * < l.-A.-.K ) Jf-s«-2 ) 

V Ai'i 2 = AiJ.'lt o ^ si i<(PH I ^ o— 1\ i o I .uPt i i I The, i\ ( j ) ) 1 
V'h'^ j = — p 1 u>i’.R / ( 1 • - X ■'^IK 1 ■*■' ( LO o/' i ^A.0 oP' I i + 
l^lAXI«slKHhi!^sI.M( THL.n{ J) ) ) 

PsLhi3L = yA|U->!' ( (/AR2 + VAR3 ) 

3C0 PRAD( K. ) =:PKAD( \ ) +PPL Auf 

2 oo Cvj.'vi r I AwS 

I t 1 1 V. vJ I’iP L L 1 U o LvJ P I 1 1 0 p\ o J I n L. i 1 I'-i L »J L R J L O f u IN 

Oi\L i-JAL.' I /-s [ 1 OM xvNvjLE 

rtRIlE{D»2i) (I»PpAi>(I)»I=l» i>i>jA<<iP 1 
2 1 F OR.'-l A T ( 4 I -jX » 1 1 3 9 3 X » 1 b 1 5 . 3 ) 1 

UbAr'iP = NbAF*P 

C F I >HU Trie I'lhAN VALUE 

P u u*''t = u o 0 

DO 400 I=i»NbA.vp 
400 PSUM=PSUM+PRAL ( i 1 
P •■iEA,\=PSUi-i/USAMP 

x.K11l(6»3u) Py.EAN 

30 r jKnn I ( / / / 2 9M .''iLAM KML)lAlt.J PKeobuiVt lL2w«b//Y) 

DO 300 I=i»NbAHP 
boo PKA I I ) =PRADU ) /Pi'i'^Ari 

!F(IXOPT,EO, 1 ) GO TO 701 
IF( N0PT,E0.4) go to 7ol 


uO TO 702 
701 CuNTIi’mUE 
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C v-ONi-'iTiON Tml 0/\TA HOR ^hl. ] PLOT 

C MAXIMUM MiMiyui'vi VALUlo IR Inc. /^RRAY 

P-iAX = -999 i Pi';iN-f?99 

DvJ 6Uo I si jiNiOAMP 

il- ( PRmI 1 ) •LT.PMlM ) PiViN = PRAl 1 ) 

1 h { PRA ( I ) oOT «P j^iAX ) Pi-'.AX = Pka( I ) 

600 COiM T I .Mot 

C plot Tht WURMnL 1 X ^ U ACOtioTlL F'pLiiouKLo 

OrtLL IPLOl I'T'Km 9 1'joAMP » L> T 9 P> i i » Ht-VvX » MiK ) 

COMTIMUL 

C PKO^lOE AKKAY C 

00 40 I=1»N5AmP 
I I=X« I 

hO UATo( I i„l)=pRADt 1 I 
MNM = 2^«MoAi'.P 

.^KI IE 16941 ) iDAlCt I ) > 1 = 1 
OaLL t'OUi< j. ( uAl CsOSAiYP s-1 ) 

Ou /O 1 = 1 >i>;6Ai-.P 

1 I I 

7 0 c( I ) =cmplx(dato( I r- 1 ) 9i; A To 1 1 1 ) ) 

U rtHRf-vY C MO* COMImIMj 7i‘L J-pLOotRCV v-OMPOMLmIo 

C i-OMPUTL rtMP L I T o lxL aMl^ Pn/iijL. '■vM'^Lt I ORiV 

»vi\irt_{694l) lCtlMl = l95o) 

FOKiiAl {iitj-oX»iEl2o39l'X»lc.il3»j) ) 

NF = M6A>‘IP/X 
00 60 1=1 jNF 
A'-iP( I }=CAoo(CC I 1 ) 

A.'IPI I }=AiV,P( I ) /Mr- 

50 A.PLI I )=2C-»ALOGlw(Ay.P( I )/2u.*L-06J 
C uji- t-rtUT .rCURiLR TR«M6rOK'-i To uAPKi-oo iN Ti-KMo or t- R c Ou:_N v. V 


C Cvj.'iPuTL rut hkht-UF.Mu Y AKKmY 

F/iAXsi./ ) 

DF = r.'iAX«2. /FLOAT t NSAi-.P ) 

Do 60 1=1 »MF 

oi = i 

60 FkFO I I } =H t*0 1 
i'(RlTLt6»6i) 

ol FORi tAT t //^4M KAoIaTLO ^aOuMu .aPtOTKuM //) 
rtPI I F ( 6»6i J t FREG t I ) soPL ( I ) 9 1 = 1 9NK ) 

62 F0RyAT(4(iFlU.2i5X9lFlu,^9bX) ) 


RH>RODtjciBTr.rry of ttt- 

OEIGINAL PAGE IS POOE 
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IFlNOPT =tQ.2) GO TO 
IF(NUPT.E^.4) GO TO 1^3 

■ jO lU IOijU 

/ijj Co,\TIi\uE 

C NOW PLOl IHc. oPli-TR ^‘1 Of" iii»- RmuImIqu jOoNo 

C ^OMi-'iTlON Inc. OATA r OR I'li. ^PlU RhL PLOl 


NF^NoAiXP/2 

CALL FPLQI IoPLsPRFQ’NF - 1 

1 UOO CON I I I'jUt 
oIuP 
L-ND 

jOtiP'-'L/ I I Nc. FOoRl t AA i A s I o I Gi'l ) 
C A ol;V|PLt. r AbT (-OoRIlIn T jI" 


0 1 1 ’l u N 0 I 0 ATA{ 2 'JAf 3 ) 

P I = AC Go I " j. « ) 

TPI=2.*PI 
1 P0“^! 
i P i = I P vj * 

I 3RLV=1 

Ou bO I3=i»IP3»IH0 
IF( I 3 - 13 PlV) 1092 o,?u 
10 TEi''Ai' = DArA(I3) 
rc.i'iPI=L)fMrt( 13+1 ) 
uAr,\ (13) =0AT Al I 3i\fV ) 

OATa( I3 + i }=DATA( 1 3RL v + i ) 

0 ATa( I 3 RLV ) =TLPPK 
0ATA( I3PL\/ + 1 ) =TL; P I 
20 IPl=IP3/2 
30 IF( I3REV-IP1) 

40 I 3 RL I 3 klV- i P 1 
IPi-- iPl /2 

IF { IPi-iPO ) i:o»3o»TL 
t)0 I 3 ,<lV=I 3RLV + IP1 
IPl-TPu 

oO IFIIP1-IP3) / 0 »lnn»lCo 
70 iP2=lPL-f:-2 

lhlLlA= rPl/FLOAl ( ISIG.\-»lP^/lPw) 
o I i>i f n = o 1 14 ( 1 h£ I A / 2 a ) 

« o ( P it = - <1 * ^ o 1 I M o I N I n 
Ao T P 1 = o I i\ t TmL f A ) 

'rt R = i « 

A I “ 0 s 

DO 90 11 = 1? IPlsIPn 
DO 80 I3=ll»IP3slP2 
I2A=I3 
I2B=I2A+1P1 

fEi'lPH = vKK'*^‘JATA(I2tJ)-Ai*L>AlA(l2L; + i) 
Tl.i'iPI=aR*uATA( I2j + 1 ) + 4l*L.»MrrtC liio) 
DATA( 12ti)=DATA( I2 A)-TEiSPK 
urt I A ( 1 Za+^ J = UAT A ( I 2 A + 1 ) — 1 Li'iP i 
DATA( I2A) =DATA( I2A )+TFMPK 
80 DATa( I2A + X ) =l)ATA( I2A+1)+1L.-1P1 


n. n 
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TEi'lPR=WR 

K = i(VKi*-W»i FPr<-/J I TP I +ivi< 

90 * i = iwl ^i'WbT J^K + T tMPK’S'/JiiT p 1 +w I 
lPi=IP2 
GO TO 60 
lUO RbTuRISi 
UiMD 

-jOBROu riNb TPLOTtpRAD? inomi'iP » Ol » PM I N * P.'lAX » i\H Ao ) 
Dli-iEMSi J,^l PKAD(llOO) s7 

c iiMi riALiz^ tml plot 

CALL PKNTON 

PLOT SIZE TO dE b «3Y 7 I-^CriL^ 
i>t\ OBJECT SPACE 

C MPmo ij lilt iNju:iULrj 01" dL'''i^<- PmooM'ji.o IC ui- PLCl I 

OPAo = .'ilpHo 
I'.-iAX = UPAO»^2u,«L)T 

CAi-L ST:32 oB( 2. s9. »2- >7. ) 

CALL wT ovJG J {!-',» I M » ktmAa ) 

C CJi'JolRuCi AaEj 

call i .ML/ 1 V 1 1 » 2 ) 

CALL AXLILI 
CALL .jAXLiT 
CALL SAXLiK 

C LABEL ^CALEG 

CALu oTiM0iv(l»2) 

CALL STNO^iCtu) 

CALL TmOoLIB 
CALL ROoLiL 

c Title iCALps 

CALL STNCn.R(U) 

C’'LL lITLi-utllH TIML (SLO) ) 

CAl,L jTfMCrtK(2yi 

^.'■vLL TiTLt-L(2'::>H l"AR ElLLi^ /vLuJolIi.. Pku.oJuKu ) 

C otNLKATE The TIMt ARTMAY 

(MT = fMPAS»2U 
T 1 =U*u 

DO 10 I = l*i\T 
^1 = 1 

10 T [ n = Tl+Dl«(UI-l,i 

c PLOT THE UATA U;ilMG ^^TKAlOiiT L h\E 6 


CALL bTNPlo(Nf) 

CALL i.LLIi-1 (TiPI^AD) 
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C TLKKJNATL PLOTTIMG 

^rtLL lXITPL 

RETuRiNj 

tiNO 

oudAOoT FFLOI ( SPL »FkL^>NF J 
D I l•lLl’Jb I uN ^PL { 1 w2a ] » FJ^E'Ji’ ( 1 ) 

C lNlllALl4c. Tiil plot 

Call ol'CCUNCAdrl 1 Q^^lL A 1 1 ui\ oPLCT^uhi 

C PLJI bI2E TO ijE 5 [3Y V I.'lCnL^ 

C JL [ uljJlC I oPaCl 

call 

CAL^ bTooO'J(Tu« } 

C CjNoTRUCT AXEb 

CALL lI,'J01V(34j»2'-'5 
CALL AXlGLI 
CALL oAXLo I 
call oAXLIk 

C .-AljLL scales 

CALL bTMDIV(3»lu) 

CALL oTMUcC(o) 

CALL XGGLiL 
CALL J'^10LGD 

C I I i Lt -iC A i-Lj 

CALL SrNCMP(15} 

CAt-L T 1 TLLti ( li:tl hRC-LENCV (hl) ) 

CALL ^lNCnK(37) 

call IIlL*-Lt3/H XFOi^TLi-J ^PL i^t. 2 ^ i * ^ 2 ) ) 

C Pi_ij| r [‘(L t-^ATA Jbij\|G oTKAiunI t-ii'tt-J 

CALL oTiNPtcj(NF) 

CALL bLLoLI (f-KEu»SPL ) 

C lEK.-iINATL PLOITING 

call lXITPL 
KL T uKI\ 

END 
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5.0 PRQGRAH SDATA 

Program SDATA is used to perform correlation and spectral analysis 
on discrete time series. The particular operations carried out on two 
given time series are: 

1) Autocovariance function for each sample and cross covariance 
between the two samples 

2 ) Autocorrelation function for each sample and cross correla- 
tion between the two samples 

3) Autospectrum for each sample 

4) Magnitude and phase representation of cross spectrum 
between the two samples 

5) Squared coherency function 

The analysis used as the basis for the program is that presented in 

SPECTRAL ANALYSIS Ai'ID ITS APPLICATIONS, by G.M. Jenkins and D.G. Watts. The 

method used is exactly that recommended in this reference. 

Two input possibilities exist. The first is to input two time series 
of arbitrary length. For this case the program will compute all the above 
covariances, correlations, and spectral quantities. The second possibility 
is to input two covariance functions. In this case the program will only 
carry out the spectral analysis. Some electronic instruments have become 
available recently which quickly compute correlation functions. This 
program can be used in conjunction with these instruments to obtain Che 
spectral representation of the data. 

A classical Fourier transform was used in the program. F’or this reason 
long samples, in terms of correlation lags, will require substantial computer 
time. Therefore, a Fast Fourier Transform program should be used for large 


samples . 
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5.1 Outlir,te of Computations 

This outline is substantially the same as that presented by Jenkins 
and Watts (pages 382 and 383) , 

Ttjo time series of data, and sampled at Increments of H sec for 

a total of K points. For convenience H is assumed equal to one. The 
frequency range is then 0 ^ f ^ 1/2 Hz. For non-unity values of H the 
frequency range is 0 ^ f 1/2H. The actual spectral levels are also to be 
multiplied by H to obtain quantitatively correct values. 

The computations are; 

(1) For the data: 

(a) the autocovariance estimate is 

i 0 i k t M-l 

t=l 

where M is the number of time lags in the estimate and the 
mean value is 

’^l W ^ ^It 
t=l 

(b) the smoothed spectral estimate using a Tukey windcw 

M-1 

= 2 (c^^(O) + 2 I C^^ik) w (k) cos -fL} 0 1 j 1 F 
k— 1 

F is a measure of the spacing of the frequency points. The 
smoothed spectral estimates are to be computed at 0,1/2F, .,.,1/2 
where F is of the order 2 to 3 times M. The Tukey window w(k) 


is defined by 
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w(k) = ~ (1 + cos “) 


|kl < M 


k > M 


(2) For the data: 

tz 


(a) the autocovariance function estimate 


'22 


1 — 
Ck) = 5- I 

t=x 


0 < k < M-1 


\jith the mean value 
N 


^2 " N J. ^2t ‘ 
t=± 


(b) the smoothed spectral estimate 


M-1 


CiiCj) = 2 {c 5 ^( 0 ) + 2 ^ c,j,Ck) w (k) cos 0 1 J 5. 


'22 


k=l 


'22 


where the same window is used (Tukey) for w(k). 


(3) For the and the x ^2 data: 


(a) the cross covariance estimate 


N-k 


=12 H ^^lt“’'l^^^2t+k ^2^ 0^k;^M-l 


t=l 


N-k 


c^iCk) - =12^-^) " N . ^^It+k 0 k £ H 1 


'12 


t=l 


(b) the even and odd cross covariance estimates 


S-llCk) - 2 '^12^“^^^ 


0 < k < M-1 


q,„(k) =y{c,„Ck) -c „(-k)} 




A f. 

' i A 


•-■j - 

/•i ":• 


l-'i 


: 

'At < 




0 < k < M-1 
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(c) the smoothed co- and quadrature spectral estimates 

2 + Z I w Ck) cos 

Q12O) ■= ^ I 

k* 1 

note that ~ ° 

(d) the smoothed cross amplitude spectral estimate 


A^^Cj) = + “^12^^^ 


0 _< j F 

1 £ j ^ F-1 


0 <_ j £ F 


(e) the smoothed phase spectral estimate 


Fi^Cj) = arctan {- = 


Ql2Cj) 


L12O) 




(f) the smoothed squared coherency spectral estimate 

Af Cj) 

q,o) " 


^iiCj) C22CJ) 


0 < j < F 


0 ±j <Y 


5,2 Correction for Bias Error ReducLlon 

T-Then there is a large delay In the correlation a large error can result 
in the coherency and phase spectrum. This bias error can be reduced by 
shifting the cross covariance function so that the maximum correlation occurs 
at zero time lap. 

The way this correction is made In he present program is to search for 
the time lag corresponding to the largest a 'lute value of covariance* A 
new covariance function is then generated such t. t the peak occurs at the 
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zero time lag. This function is Fourier transformed and the phase shift 
Introduced by the time shift is applied to the spectral estimate to correct 
it back, A consequence of the generation of a new covariance function is 
that it necessarily will contain fewer terms than the original. This results 
in a wider frequency bandwidth in the spectral estimate. This is why the 
program will seem to generate different bandwidths for different batches of 
data even when the same number of lags are called oiit in the input. 

5.3 Input to Program 

Two possii'-ilities exist for input to this program. The first is to 
read in two time series of variables for which correlation and spectral 
computations are required. This is obtained from data card one by setting 
the variable OPTION = 1, FORliftTClIS) . If OPTION is not set equal to one the 
program will read in the data as if it were covariance functions. Input 
for OPTION = 1 is described in Section 5.3.1 and for OPTION 1 in Section 
5.3.2. 

CARP 1, FORMAT (115) 

OPTION = 1 read in time series 

^ 1 read in covariance functions 

5.3.1 OPTION = 1 
CARD 2, FORMAT (215) 

NSAMP number of points in the time sample 

MLAG number of time lags desired in the covariance functions 

typically MLAG is chosen to be approximatelv one tenth as large as NSAMP. 

CARD 3, FORftAT(lFl0.2) 


DELT 


time increment in samples 


(Sec) 
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The remaining Inputs under this option are the actual values of the variable 
in the time series. They are read 7 values to each card in the FORMAT (7F10. 2) . 
The number of cards required depends upon the size of the sample. The com- 
puter call is 

DO 10 I = 1,2 

10 READ(5,2) CX(J,i:,I = 1,KSAI^) 

Tills is the total input required for this option. 


5.3.2 OPTION 4 1 
OARD 2, format (2F10, 2) 

XLAO number of values in the covariance function samples 

DULT time increment in the functions (See) 

The remaining cards are used to read in the covariance functions. Tn the 
present format the functions are read In as follows. 


DO 10 I ^ 1,2 
DO 10 J = 1,2 

10 READ(5,2) (COV(K,I,J), v ^ i^Mi) 


where Ml = XLAO + 1 and the 
C0V(K,1,1) = = 

COV(K,2,2) = C^jO') = 

cova:,i,2) = 0^2 00 = 


indices I,J mean the following: 

auto covariance function for series 1, k = 0,M 

autocovariance function for series 2 , k = 0 ,M 

cross covariance function for positive time 
lags, k = n,M 


C0V(K,2,1) = 0 - no = cross covariance function for negative time 

lags, k = 0,M 


The format used to read in these functions is FORMAT {7F10. 2) . This completes 


the input requirei^nts for this option. 


5,4 Program Input 


Tlie output of the present form of this program is simply a listing of 
each of the functions computed. The list of functions computed was pre- 
sented in Section 1.0. 


nnn n r r n nnnonr. nn 


5.5 Program Llstinp; 
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PkOGRA.’^ 5DATA( InPUT »UUTHvj1 , I AHE3=INPoT s T ''.PE6 = ^'0 I HuT ) 

M STmTISTIv^AL UhTA AfSjMLYiii:^ PKOuRnM wrliCM COi'lF^vJ I Lo F ni_ rOfSlt-iiONi> 

AUTO COVhi'J i ANlIi AfyjU GROij'i CvjVr\j<Jl MNGc A\S «Jc.LL 
AOTO hNl> -iPFCtKuf/i 

RtAD IN the tike SERIEo 

NbA.'iP = N^KBEK OF TI-'.E b^KPLEb 

KLAO = NvJMulR Ol' CORRELAllON Lwbb 

OELT = Tl.'iE INCRcvEnT (bLCJ 

AlI»J) = I HE I vvO SAi'tHLE-J lo DC. MNAi_Vi;.EO 

DATA PI/3.14159/ 

READ(5»1) OPTION 

option = ^ RtAu TJ-KL bcRiLb QfHiiRk'VloL RlmlJ COv^RiANOLb 

IF (OP T I ON 0 to. 1) GO TO 5u 
READ(5»i) NbAi'iP »!mLAG 
KtAO(5»2) DELT 
DO 10 I=i»2 

10 READ(5»2) ( XU» I 1 » J=1 »NoAHP) 

wRi TE( 6»40 > ( I »X (I » i ) »A(1 ) » 1 = 1 »NbAMP ) 

40 FORMAT ( 2( 5X»1I5»5X»2E?^.51 ) 

Form eov^-'KImN'-ls anl^ cokkll«i ^ons 

CALL CuV AR ( NoARiP 9 X » COV » inL Ac ) 

OQ,-'iPuTt Inc_ oROb;j POifttR oPLClR^L Ui..NoiTY 


GO TO 22 

50 READ! 5»2) XLAGjO'-LT 
mLAG=IFIX(XLAG) 

M1=MLAG+1 
DO 2i J=i»2 
DO 21 I=l»2 

21 READ (5 >2) ( COV ( K » I » J ) » K= 1 9 K 1) 

22 CONTINUE 

COVaR I ANC L.:^ ARE RfAl^ IN T nb ORUtR t 1 * 1 1 » ( 1 » 2 ) > ( 2 » 1 ) » I 2 ’ 2 ) 

'wHLL CRbPEC (COV»MlAu »ucL 1 »N>:'^iMP ) 

1 FURHATC2I5) 

2 FORi^iAT ( 7 F 1^«2 ) 
oTOP 

END 

^juBRUuT I NE CUVAR t NbA'’tP 9X1 CU\/ »i’ilAG ) 

UlfinNblQN bM(2) »X ( 1 J24» 2 ) 9CuV I iui » 2 » 2 ) 9 uCO V t 1 u 1 9 2 9 2 ) 9 
lOJRt 1 u 1»2 92 ) jOCIja ( 1^1 92 92 ) 9 bur. ( 2 ) »xi'. ( 2 ) 

The muTU ^ND CRO.aS CvJVaK i ANLtio AC.LL /\o Tnii auIO AND 
CKObb correlations ARE CuKPuTeD 
COMPUTE MEAN VALUES 


EEPEODUCBJLnY OP TIP 
ORIGINAL PAGB IS POOR 



00 1 J=l»2 
i>UM{ J)=0o0 
DO t I=1*N^AMP 

2 J)t=SuiVi( J)+X( I »J) 

1 Xi'i ( J ) =bJ*''i IJ ) /Ni>A<*iP 

I 

C PRIMT MEAN VALUED 

rtRITE(6»10) 

10 FORtMAT ( iH-L * JbH RES'JLiij OF COKj^ELaIIuN Cui'iPu UT I oNo 
A'RlTEt6»ll) (XM( J) ♦ J=1 »2) 

11 FURKAT ( //^14H MEAN VALUE oEKlto ONE = iE^w.b / 

12An MEAN \/ALOE btRluo T ViO = ic.2'-'»b //) 

C COMPUTE ui. VI AT IONS r ROM iHd MLrtN 

DO 3 J=l»2 
DO 3 i =5 Is NS AMP 

3 A ( I » J } =X ( I * J ) -Xm t J ) 

C CUf'iPUTE CuVAKIANCF^j 

M1=MLAG 
DO 4 X= Is i'll 
DO 4 J=ls2 
DO 4 L=ls2 

b U iM 1 = u B U 

('.d = NbAMP— s 
uO 3 I-lsi-i2 

5 ^J.'il = X(IsJ)*X(l+N-lsL )+bUl«';l 

4 COV( As JsL)=bUMl/FLOAl (NbA.^, Pi 

I ’ C COM.PUTE DIFFERENrFD COVmRIANCEo 

I 

j DO 6 L=l»2 

I DO 6 J=ls2 

OCuV { 1 s Js t ) =-C0V (2*‘-sJ)4-c.'*^Cuv{l>Jsi-)-CuV(2susL) 

^ <vt3 = Ml-l 

I DO 6 K=2si''i3 

^ 6 DCUVI JsL ) =-cOv I K-1 » JsL) +2«*Cuv ( As JsL ) -CU\/ (A + 1 s J sL ) 

E WRITE{6»12) 

1 12 FORMAT! 3in TABLE oF AUTO COVARIANCE <1-1) ) 

■ vuKIl£(6sl4) (X »COV ( K.S 1 s 1 1 sN=i s,ni ) 

j 14 FORi'iAT (4 ( IIS^'sXsIpiS.SjSa ) ) 

rtRlTE(6slb) 

! 13 FORMaT(J1*i TAdLE oF muTO CO'/^RIaNCL (2-2) ) 

i aRiIe(6»14) I a sCOV C n s 2 »2 ) »^=l »‘Al } 

i rtRITE(6sl6) 

i 16 FORMAT (32M TABLE OF CROoo COVARlANCt (1-2) ) 

4 WRlTE(6sl4) ( K sCOV C X s 1 s 2 ) s^=l sMl ) 

i . rtKirE(6sl7) 

1 17 Format (32i1 table of CKUSb CoVAkIAivCE ^d-l) ] 

^ aRITE(6s14) (XsC0V(As 2»1> *A=l»'’'a) 


OUTPUT DlhFERENCED COVARlANCtib 


o ri n n ri r. r. r. 
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wRITE{6»181 

i8 hORMMT ( itix j46h Ui FE cRh.N'-L.u X.OVmK loAMi- ORi^tiR Ao AdOVc.) //) 

CDHPUTL LOkRuLAT lONii NORi'IaL A COVmR i AN(-lii> 
ixOTn normal AND u 1 F t* cRtN^cu rORiviS 

DO 20 L=l*2 
DO 20 J=l»2 

COVn=bQRTlCOV{l»J,J)*cOVtl»L.»L] ) 

DCOVMsSQRl (DCGVd ,J»J)*DCOV(1 »L»l'i 5 
DO 20 K=1»M1 

C0RtK9j»L)=C0v^ J » L) /COvM 

20 DCOR{ J»L)=DCOV( K» J»L}/l.COVM 

PRINT Thl LORRLLatIONS ANl, IriiiiR u i |- h lRuNLuU t*OKiMS 
wRI TE(6»21 ) 

21 Format ( Itii »50n auto ani^ v.rq:30 uorKl-laF iCNS anu l/ i n- c.Rt_NUiii^ rOHMS 

i///4ut-i AUtO I-ORRI-lAT iON diln U i, F r lRuNUlu rORM) tl-1} //) 

rtRlIElOslA) tKsCoRd^l*!) 
aKITE{6922) 

22 FORMA! { //£;3h AUTu CukRtLMTlU(\ ) 

aRI rEl6slA) 1?;9 CoRU»2»2) ) 

wRITEt6,23) 

23 Format ( /// 25ii cross corkl-latiun (i-2) ) 

rtRJ TE C6f lA) (RjCOR ( Nd >2 ) ».'v=i ».-a ) 
rtRlTLt6i2A) 

24 Forma 1 ( / / / 2 5ri CROSS CuRKlLA I loN i£.-l) ) 

«(iRlTLt6»l4l lrs9CCR^R92»l)»A-l»'''il) 

Trld COKPLcTlo TME CORRELA f Ji Oivi PfIaol 

RETURN 

END 

^JCjRUUI iNu CR^PLC ( CUV 9MLAU»ut-Ll » NohMP ) 
r 1 iE CK^ jo Fu*lK jpEClKAU Ot.i'i'jilt lo Coi’ih'Ulc.L' 

D 1 1*. EK.J i '-'N CuV(Ivj 152»A) sC'^VPli'ji) »oPtCP!3^3) *oPl.l( 3‘-'3*2) sFRE'^^tSuS} > 
lEVtiOll »X'^Dtlul)»SO(ju3j9 CuoP lC 13^3) »'2(oFEC(3k^31 »i^i-tAs:»E(3'^3) » 

2CuH:30{ 303} jASPEC: 303 »2) »uuMl^u2) ,UUPt2o2) 

DIMENSION CSPEC(3o31 »FKEC^3^^J) 
data PI/3. 14139/ 
i''ii = i'‘iL AG 


F i b rtiiuRc. b is TltE N^iMixi-ri OF liMt- LAub 10 nLiUN f fit- 
FROx.t.jbLo ft'ITFI .'lAXIi'iuM CoRRELA ! i Oi'l AT 2ErvO 
UETbRMlNb LAG TIME FOR ScRUb alIoN-IcNT 

rtRI TE ( 6 *7U ) (CUV ( I » 1 »2 ) »CuV 1 1 i2 d ) d =1 *M1 ) 

DO 21 K=1»MlAG 

21 DbMlA)=CUV(Ml + l“|ivj2d) 

DO 22 R=2»isl 

22 DUi-lt i’tLAG+N-1) "COV ( K.d »2 ) 

MM=2*MLAG-1 

•vRi TEt6»26) ( I »DUM( I ) d = l 


ill':FR(A')b\3oVj'[.l i’Y OP THE 
ORTOO.AL rs POOR 


n r 
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26 F0RMATt4OX»lI5»5X» ) 

XMAX = -999<. 

00 2 3 K = l»iSM 
0JP(O=ABotDU,'^(K) ) 
ir(DUP( N) .Gr.XHAX ) GO 24 
GO TO 23 

24 XMAX = DUP(i<) 

IG=K-M1 

23 CONTINUE 

iVRITE(6j2P) IS 

25 FOHi'iAT ( ///22H ALIgNi'.CNT LAG -W.-iuEN lllu///) 
^ONultlON 1 UL UM I A I- OK Ifii- OrtLCuL^'iiON Of Thl 

CKOob GPcClKi^i'i 

M2=Kl-IAiib( Ib) 

^ONL^^TiON Tiiii OrtTA fOR ALrTObp;-^-T Ro,X LrtLLuLf->T iON 

DO 20 K. = i9C 
00 10 I=l>i'il 

10 COVPC I I ) 

CALL A^JopcC ( I'ILAu »:iELi 9 CovP»»jHl.CP»i-Klws lo) 

C oTJKc AJfJbPtCTRUM 

JN F = 2 ■5'" i'I2 + 1 
00 11 I=i>NF 

11 oPLC( I»rv)=oPLCPl i ) 

2 0 CONI I No 1£ 

c.V( 1 ) =Du;'WM1 + I3 ) 

XOD ( 1 ) =0 ,0 
DO 40 K = 2»»X2 
rl ''i = i'^ 1+ I o + N — 1 
Hi'l(‘iri = i*(l + I b — X,+ 1 

E\/tK)=Oo 5 '^ ( OOi'i I .‘I *iM ) +L/ 0 .'-i i ) ) 

40 X uO ( K ) = 0 ■ 3 ^ ( OOj'. [ >‘1 ) ~ Do i‘i ( I ) ) 
fiKirn.{6s7^^J lE''ll)>AODll)»I = l»''i2) 

c 00 The fo^rilr ikanokoKi'S 

CALL CRCbPEC < LV »A0D*C0bPt,C 9vjbHLL ».-.LAG »0LLT » i o ) 

C Cu.'lPOTE PhAGL ANolL ANo CoMLklNOE FuNCI IuN 


v_ Kt-IijKN ^KObG bPLLyRvjV. To ONofiiFItu fOKi'i 
DO 60 I=lsNF 

AR>j = P1*Fl0aT ( i-1 )4tFL0AT ( ib)/FL0«r(NF ) 
CN=CObtARG) 

■^N = J I N ( ARo ) 

oU|Vll = CObPi-CC I )#bN+GbPEC( i 
OoN2 = C0oPcC( i )*CN-Q6PEC( 1 
COSPECt I ) =DUN2 
OSPECI I )=0UM1 
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iQ(I J =COiJptC ( 1 ) *CoSPt.C( I ) +Qi>PLL ( I )-HQi,pLC( 1 ) 

IFCCOSPEC( I) .EQ«0*0) Gu TO bU 
PHAbt [ I ) =ATAfvl2 C-oiSPECJ I ) »C0oPEC( I ) ) 

GO ro 81 

80 PHAGEt I )=Pl/2. 

81 CONTINUE 

PHA^EC I )=PHAGE( D^abo./pl 
CGPEC( I )=GQRT(SQ{ T } J 

CUHoOd )=ou( I)/(oPEC( I»l)^^oPEC(I »2) ) 
if(gpecci»1)«le.u,0) go 10 51 
GO TO 55 

52 CuNTlNuE 

IF(s>PEC( i »2) .lE.^«u) go to 53 
GO To 54 

51 ASPhC( I .1)=-100. 

GO TO 52 

53 AoP^Ct I *2)=-100. 

GO t xJ 60 

5 5 CON T I NUE 

AoPECn .1 )^ALOGl^(SPEC( I >1 ) ) 

GO ro 52 

54 CONTINUE 

AoPEC{ I »2 ) =ALuGlu (SPEC( ! *2) ) 

60 C'Ji'i T I Nu£ 

C uu Ni-RATE The FKLOUENCY AKKAY 

•jl=±,33/ (FLOAI (i‘.LAG)*DEL1 ) 

UO 4i I =1 »NF 

ri ( Rl.u( I )=FL0AT( I-l)/(2««'L;LLT'*fFLiJ<AT(Nf'-l} ) 

O rtRlli- THc_ ^PLCTKAl LLVfiLO (^-KO^o PcaO) jPhAoL mNu L(jl(tKLK*^Y 

/IRITE(6»61) 

61 FcRMmTCIhi »35fi Rc-SULlo 01“ PO»lR oPcl I RAl AN^ltoIo /// 

IIVH AoTOoPECfRuH (1-1) //) 

70 FuR.’iAK 3(jA9lFl^«4»5XsiEl5.G)) 

rtRI Tu(6»7o) (FRLQ( I) » GPl C ( I d ) » I = 1 »NF ) 
tvRl TE(6»6<^) 

62 FORHAT ( ///19H A0TCGPEc1R'^>'i 12:“2) ///) 

aRI TE (6 *7u) (FREq( n»SPhC(I»2)9l=l»NF) 
aRITE(6»65) 

65 FORHA ! ( ///15H CROGg oPECRuP ///] 

»vKI TE ( 6*70 ) ( FKEQ ( I ) *C5HEC ( I ) » I =1 »iMF ) 

rtRITE(6*63) 

63 FoKmATC /// 12H PHASE ANQLt ///) 

«iRI IE(6»7U) '(FRtQ( I ) *PHMot( i } »I-i»NF) 
rtRI TE(6*64) 

64 FORiMT ( ///IIH coherency ///) 

rtRI TE (6*70 ) (FHEQ( 1 ) >C0HSQ( I) » 1=1 *NF) 

return 

END 

ouBRiJoT INt CKuoPEctEV * XOu 9COjPto*woHLC»nLAG»DEL] »Ii) 
l) I r^tNS I ON EV(lUl) ,X0u(lL.lJ »Cv^oPlC I 3v^3 ) »OoPEC(3u:j}*rt(lul) 

L COMPUTES OMOOTHEO DO- ANO QUmU- SPclCTRm aNO oOUf^RLuMMPL 1 f U ul 


n n 
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DATA PI/3,1A159/ 
i^il = MLAG 

NF = 2^tM+i 

C COMPUTE wuIGHTS L»Sli’'lU TunEY vVlNuOk^ 

DO 20 J=1»M2 

2 0 w( J)=0. 5* ( 1, + COotpi^^FLUAl t J ) / FLOAT I i-i) ) ) 

DO iu I=1»MF 
O U 1*1 1 = 0 e 0 
ooKi2 = 0 « J 
DO 30 K = i»yi2 

AKG = PI^tFLOAT lf;)*FLOAT I I-l (/FLOAT (NF-l ) 
oN = 0U4CAR0) 

CG=COS( ARG) 

o»jt’ii = iUMl + E V I fv+ 1 ) ■‘■■(V t fs ) *CO 
30 oo,n2 = boiM2+XOD( f^+1 ) ^ 

COoPtCll )=2.*DEU*{tL\/m+2.*ouivil) 

10 ObPcC ( 1 )==4.ttDELT*SUM^? 

KLTuRN 

Ei^jD 

oOBK^o T I NE AO op EC ( ML AG » DLL I » Cov »dPEC»FKL*J* I->) 

C COMP'jTL AuTOiPECl rum Oh CQVmKImKClo 

Dl.-lLNolUN COV (lvl)j*(lolloPtCt3J5) » FRlO 13^3) 
DATA PI/3.14139/ 

The tukcY w1i\doa is used fur data c5i''iuurHiNo 

CuMPUr WEIGHTS 

Mi-HLAG 
iM = .Ml-IABul IS) 

,\[F=2*M+1 
DU iO I=1»M 

10 M I I )=u.3*t lo + COo(pI*FLOAl U )/FLUAltri ) ) ) 

C CrtLLoLATE H smoothed AuTOSPLCIRml lsIIMATE 

MM=M-1 

UU 20 1=1 *NF 
o U Fi = 0 . u 
DO ^11 K = l*iVi.M 

ARG = PI^tFLOATI K)*Fl 0AT t I-l ) /FLOAT ( NF-1 ) 

VAR = C0V(K:+1 )^^W(K)*COS(aRG1 
21 sUM=SUM+VAR 

20 SPECt I 1 =2.*DELT^lC0V( 1 )+2,*SUiv.) 

C COMPUTE D'^NDwIdTH ANU DhuRELo oF hREEuOM 

return 

END 


